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ABSTRACT 



i Breakthrough direct detections of planetary companions orbiting A-type stars confirm the existence of massive planets at relatively 

large separations, but dedicated surveys are required to estimate the frequency of similar planetary systems. To measure the first 
estimation of the giant exoplanetary systems frequency at large orbital separation around A-stars, we have conducted a deep-imaging 
^ C"| survey of young (8^100 Myr), nearby (19-84 pc) A- and F-stars to search for substellar companions in the ~ 10-300 AU range. 

O { The sample of 42 stars combines all A-stars observed in previous AO planet search surveys reported in the literature with new AO 

observations from VLT/NaCo and Gemini/NIRI. It represents an initial subset of the International Deep Planet Survey (IDPS) sample 
of stars covering M- to B-stars. The data were obtained with diffraction-limited observations in H- and A" s -band combined with angular 
differential imaging to suppress the speckle noise of the central stars, resulting in typical 5<x detection limits in magnitude difference 
of 12 mag at 1", 14 mag at 2" and 16 mag at 5" which is sufficient to detect massive planets. A detailed statistical analysis of the 
survey results is performed using Monte Carlo simulations. Considering the planet detections, we estimate the fraction of A-stars 
having at least one massive planet (3-14 M Jllp ) in the range 5-320 AU to be inside 5.9-18.8% at 68% confidence, assuming a flat 
distribution for the mass of the planets. By comparison, the brown dwarf (15-75 Mj up ) frequency for the sample is 2.0-8.9% at 68% 
confidence in the range 5-320 AU. Assuming power law distributions for the mass and semimajor axis of the planet population, the 
qq , AO data are consistent with a declining number of massive planets with increasing orbital radius which is distinct from the rising slope 

. inferred from radial velocity (RV) surveys around evolved A-stars and suggests that the peak of the massive planet population around 

^ ' A-stars may occur at separations between the ranges probed by existing RV and AO observations. Finally, we report the discovery of 

' three new close M-star companions to HIP 104365 and HIP 42334. 

' Key words, instrumentation: adaptive optics - instrumentation: high angular resolution - methods: observational - methods: statisti- 

, cal - stars: imaging 

(N 

1. Introduction ets in the Solar System, remain outside the range of detection 

, . , , ,. , , of these methods, and, consequently, little is known about the 

^ ; An extensive population of exoplanets has been discovered down extrasolar popu i ation of wide orb it planetary systems. 
5_J to sub-Jovian masses and at orbital separations below 5 astro- 

nomical units (AU) based on large-scale radial velocity (RV) Direct detection with high-contrast imaging provides a 
surveys (Mayor & Udry 2008; Marcy et al. 2008) and on tran- method tQ explore the wide orbk planet population . Several 
sit surveys, which are now identifying hundreds of new candi- adaptiye Qptics (AQ) suryeySi concentrating on Solar-type stars, 
dates (Cabrera et al. 2009; Borucki et al. 201 1). These indirectly haye been conducted to search for i ow . mass subs tellar compan- 
detected planets provide invaluable information on the distnbu- [om around nearby young stars Becmse direct imaging needs 
tion of close orbit planets (Cumming et al. 2008) and on their tQ overcome the large contrast ratio between the star and a 
frequency around nearby stars covering a large range of masses potential subst ellar companion, the existing surveys were fo- 
(Marcy et al. 2008; Johnson et al. 2010a). Planets at orbital radii cused on FGKM stafSi around which young objects down tQ 
larger than 5 AU, compar able to the locations of the giant plan- a few masses of Jupher (Mjup) wouM be detectable at sepa . 

* Based on observations collected at the European Southern ration lar S er ^ an a few tens of AUs ' The ma J orit y of these 

Observatory, Chile, ESO programs 081.C-0519, 083.C-0706, 084.C- surveys reported no (Masciadn et al. 2005; Billeretal. 2007; 

0605, 087.C-0559, 088.C-0477, and at the Gemini North observatory, Lafreniere et al. 2007a; Kasper et al. 2007; Leconte et al. 2010; 

Gemini programs GN-2007B-Q-59, GN-2008A-Q-77, GN-2008B-Q- Janson et al. 201 1; Delorme et al. 2012) or few (Chauvin et al. 

64, GN-2009A-Q-80, GN-2009B-Q-93. 2010) substellar companions to nearby stars and have placed 
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upper limits on the population of massive planets at large 
orbital separation. Some objects were nonetheless discovered 
(e.g. Lafreniere et al. 2008; Thalmann et al. 2009; Billeretal. 
2010), demonstrating that these objects exist, but are indeed rare 
(Nielsen & Close 2010). 

Recent breakthrough detections around young A- 
stars - HR 8799 (Marois et al. 2008b, 2010b), B Pictoris 
(Lagrange et al. 2009, 2010) and Fomalhaut (Kalas et al. 2008; 
Janson et al. 2012) - have provided new perspectives on 
the search for companions orbiting more massive stars. The 
results of RV surveys of the evolved counterparts of A-stars 
(Johnson et al. 2007) have also provided intriguing results 
(Johnson et al. 2010b; Bowler et al. 2010a; Johnson et al. 201 1), 
with significant differences in the planet population compared 
to RV-detected exoplanets around lower mass stars. Statistical 
analysis of the results clearly shows a higher frequency of 
planets around more massive stars and larger masses for the 
detected planets (Johnson et al. 2010a). Moreover, theoretical 
work on planet formation by core-accretion similarly shows 
the same kind of correlation between the star and planet mass 
(Kennedy & Kenyon 2008; Alibert et al. 201 1). 

Although more technically challenging in terms of obser- 
vations, these recent results suggest that more massive stars 
may present more favorable targets in terms of exoplanet de- 
tections. And estimations of planet yield show that A-stars will 
be even more favorable for future large direct imaging surveys 
(Crepp & Johnson 201 1) with dedicated upcoming instruments, 
such as Gemini/GPI (Macintosh et al. 2008) and VLT/SPHERE 
(Beuzit et al. 2008). To pursue the possibility of a higher fre- 
quency of massive planets around early-type stars, we have con- 
structed a sample of 42 young A-F stars observed at high con- 
trast to be sensitive to massive planets. We report new measure- 
ments for 39 stars, and we include 3 A-stars from the literature 
that have been observed in previous surveys. This survey intends 
to start defining the population statistics of massive planets and 
brown dwarfs (BDs) at orbital separations in the tens of AUs 
from their parent A-stars. This sample is an initial subset of the 
International Deep Planet Survey (IDPS) spanning M- to B-stars, 
the results of which will be presented in a forthcoming publica- 
tion (Galicher et al. 2012, in preparation). 

The selection of the target sample is explained in Sect. 2. In 
Sect. 3 the observing strategy, observations and data reduction 
steps are detailed. The detection limits of the survey are derived 
in Sect. 4, and we describe the identification of the candidate 
companions detected in the data. In Sect. 5 we use the detection 
limits to perform a statistical analysis of the survey, from which 
we derive a first estimation of the planetary systems frequency at 
large orbital radii around A-stars. Finally, we discuss our results 
and conclusions in Sect. 6. 

2. Sample selection 

The sample includes all A-stars previously observed in high con- 
trast AO imaging surveys sensitive to planets (Nielsen & Close 
2010; Chauvinetal. 2010; Janson et al. 2011) combined with 
the 35 A-stars and 4 F-stars observed for this study, resulting 
in a total of 42 stars, with 38 A-stars and 4 F-stars. Of the nine 
A-stars included in previous AO surveys, six were part of our 
selection (see below); they were re-observed and the new obser- 
vations provided deeper limits. We used the existing surveys and 
the new observations to define the target set listed in Table 1, 
and used the best limit for each target. The sample A-stars span 
the full A0-A9 spectral type range. Because Fomalhaut was not 
observed as part of an AO planet search survey, Fomalhaut b 
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Fig. 1. My vs. B - V color-magnitude diagram of the A-stars of 
our target sample. The stars with published ages are plotted with 
a green circle. Stars without a published age are noted with blue 
triangles and have been assigned an age of 125 Myr, based on 
their positions within the region defined by the Pleiades A-stars, 
which are indicated with orange stars. An empirical linear fit 
to the Pleiades is also shown. The older population of 625 Myr 
Hyades A-stars is plotted with blue squares. 
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Fig. 2. Age, distance and spectral type of all the stars in our sam- 
ple. Table 1 gives the individual properties of each target. The 
size and color of the symbols correspond to the spectral type of 
each target: bluer colors for earlier spectral types and redder col- 
ors for later spectral types. The points grouped at 125 Myr cor- 
respond to stars without a published age, as described in Sect. 2. 
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(Kalas et al. 2008) was not included in the sample or the subse- 
quent statistical analysis. 

In order to define a representative sample of young A- 
stars, the A-stars newly observed for this study were selected 
to have positions on the color-magnitude diagram within or be- 
low the band associated with A-stars in the Pleiades. For the 
A-stars in the sample, the positions of the stars on the color- 
magnitude diagram are shown in Fig. 1. The target ages di- 
rectly impact the minimum detectable mass, given the mono- 
tonic decline in planet brightness with age (e.g. Fortney et al. 
2008; Chabrier et al. 2000). The majority of the A-stars - 22 tar- 
gets - have age estimates from the literature combining several 
techniques (Rhee et al. 2007; Moor et al. 201 1; Zuckerman et al. 
201 1; Janson et al. 201 1; Moor et al. 2006; Chauvin et al. 2010; 
Zuckerman et al. 2001; Tetzlaff et al. 201 1; Stauffer et al. 1995). 
These values are reported in Table 1 and used in the analysis of 
mass detection limits. Of the 22 A-stars with literature ages, 10 
are members of nearby moving groups or associations of young 
stars. For the remaining 17 A-stars, the age is based on a com- 
parison of the position of the targets and the best-fit to the dered- 
dened Pleiades A-stars, a population of ~ 1 25 Myr (Stauffer et al. 
1998) A-stars. The individual reddening values were taken from 
(Breger 1986), and Pleiades A-stars on the photometric binary 
sequence were excluded from the fit. As shown in Fig. 1, the tar- 
gets without literature ages are located generally very close to 
the Pleiades fit and within the range of Pleiades members, so an 
age of 125 Myr is adopted for these 17 stars. To ensure accu- 
rate placement on the color-magnitude diagram the A-stars were 
limited to stars with Hipparcos parallax uncertainty of <5% 
(Perryman et al. 1997). A distance cutoff was also imposed to 
probe angular separations corresponding to orbital separations 
similar to the outer regions of the Solar System and within the 
radii of protoplanetary disks imaged in scattered light and mm 
emission (e.g. Fukagawa et al. 2004; Mannings & Sargent 1997) 
around Herbig Ae stars, the precursors to massive stars. Both 
debris disk systems and stars without detected excess emission 
were included. Of the 38 A-stars in the sample, 14 have excess 
emission from debris disks, as noted in Table 1 . 

In addition to the A-stars, 4 early F-stars with F0-F5 spec- 
tral types were observed. Since F-stars do not evolve off the 
Main Sequence as rapidly as A-stars, the color-magnitude di- 
agram is not as a reliable tool for age assessments of young 
F-stars. Consequently, the F-stars were limited to stars with 
age estimates from the literature based on space motion linked 
to nearby moving groups or associations (Mooretal. 2011; 
Zuckerman et al. 2001; Rhee et al. 2007). Of the 4 F-stars, 2 are 
debris disk systems. A graphical summary of the sample popu- 
lation incorporating spectral type, age, and distance is given in 
Fig. 2. In summary, the sample consists of a total of 42 stars, with 
38 A-stars and 4 F-stars with a median age of 100 Myr and a me- 
dian distance of 50 pc; key parameters of the sample are listed 
in Table 1 . All members of the sample have no catalogued visual 
binary within 5" to ensure a detection limit that is uniform with 
position angle. The targets with wide (>5") common proper mo- 
tion companions from the WDS catalogue (Mason et al. 2001) or 
close (<0.1") speckle companions (Mason et al. 1999) are noted 
in Table 1 . 

3. Observations and data reduction 

Data were acquired between 2007 and 2012 with VLT/NaCo and 
Gemini/NIRI, two instruments located at the focus of 8 m tele- 
scopes. NaCo is installed at the ESO Very Large Telescope in 
Chile. It is the combination of the NAOS adaptive optics sys- 



tem (Rousset et al. 2003) and the CONICA observing camera 
(Lenzen et al. 2003). It provides diffraction-limited images in 
the near-infrared on a 1024 x 1024 Aladdin InSb detector array. 
Our data were acquired using the S13 camera, which provides 
a spatial sampling of ~13 mas/pixel (better than Nyquist) and a 
field of view (FoV) of 14" x 14". The targets were all observed 
with the broadband ZT S filter (2.00-2.35 pm). 

NIRI (Hodapp et al. 2003) is installed at the Gemini North 
telescope in Hawaii. It was used in combination with the 
ALTAIR AO system (Herriot et al. 2000) to obtain diffraction- 
limited images on a 1024 x 1024 Aladdin InSb detector array. 
The f/32 camera provides a sampling of ~22 mas/pixel and a 
FoV of 22" x 22". To improve the Strehl ratio and provide 
sharper image quality at large distance from the star, ALTAIR 
was used with the field lens inside the optical path. The tar- 
gets were observed either with the CH4-short filter (CH4s, 1 .54- 
1.65 pm) or the K' filter (1.95-2.30 pm). 

Both instruments were used in pupil-stabilized mode (pupil- 
tracking mode on NaCo) to reduce the number of independently 
rotating speckle patterns, improve the stability of the point- 
spread function (PSF), and provide rotation of the FoV, which 
is then used for proper implementation of angular differential 
imaging (ADI; Marois et al. 2006). To optimize the detection 
of faint companions, the science images were saturated by us- 
ing detector integration times (DIT) between 15 and 30 s in 
most cases. These DITs allow the detections to be background- 
limited at large separation, but since our targets are very bright 
(H < 7.0), the regions up to 0.4-0.8" are in general completely 
saturated and have been masked in the analysis. Each observing 
sequence started with the acquisition of an unsaturated PSF at 
five dither positions, obtained with the ND_short neutral density 
in NaCo (transmission of 1.20+0.04%) or a short DIT combined 
with narrow-band filters in NIRI (//-continuum in //-band or H2 
in ZT-band). These PSFs were then followed by 100-250 deeply 
saturated science images. For some observations, sky frames 
were acquired prior to the science frames following a dither- 
ing pattern. Due to restrictions on detector saturation, second 
epoch data on NaCo were acquired with a much smaller inte- 
gration time and a larger number of exposures. Some NaCo data 
was also acquired with a much smaller number of exposures but 
using co-addition of 5 exposures. All the observations are sum- 
marized in Table 2. 

The data reduction followed standard procedures. For the un- 
saturated dithered PSF frames, a sky frame was constructed by 
taking the median of the images. The individual images were 
then sky subtracted and divided by the flat field before being 
registered to a common center. The final unsaturated PSF was 
obtained by median combining the five images. This PSF was 
used for the photometric calibration and calculation of the de- 
tection limits. For the saturated images, a sky frame was sub- 
tracted when available, and each image was divided by the flat 
field. Bad pixels identified from the flat field were replaced by 
the median of adjacent pixels using a custom IDL routine. Field 
distortion was taken into account for NIRI since it reaches sev- 
eral pixels in the outer part of the field: all images were corrected 
for this effect using a custom IDL routine and an optimized dis- 
tortion solution that will be presented in Galicher et al. (2012, in 
preparation). For NaCo, the field distortion was found to be neg- 
ligible (< 1 mas) with the S 13 camera in the 7" x7" inner part of 
the detector (Trippe et al. 2008; Fritz et al. 2010). The center of 
each saturated image was estimated using the fitting of a Moffat 
profile (Moffat 1969) on the PSF wings. The saturated region of 
the images was assigned a zero-weight in the fit to avoid any 
bias. For NIRI, the center of the saturated images moved by up 
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Table 1. Target sample and properties. 



HIP 


HD 


HR 


a 


6 


Distance 


Sp. Type 


Age 


Age ref . 


MG/assoc. 


B - V 


V 


H 




IR excess 


IR excess ref. 


Binary 








(J2000.0) 


(J2000.0) 


(pc) 




(Myr) 






(mag) 


(mag) 


(mag) 


(mag) 








Observed targets 


7345 


9672 


451 


01 34 37.8 


-15 40 34.9 


59.4 


A1V 


20 


1 




0.07 


5.62 


5.53 


5.46 


yes 


1 




10670 


14055 


664 


02 17 18.9 


+33 50 49.9 


34.4 


A1V 


100 


1 




0.02 


4.03 


3.86 


3.96 


yes 


1 




11360 


15115 




02 26 16.2 


+06 17 33.2 


45.2 


F2 


12 


2 


P Pictoris 


0.40 


6.79 


5.86 


5.98 


yes 


1 




12413 


16754 


789 


02 39 48.0 


-42 53 30.0 


35.7 


A1V 


30 


3 


Tuc-Hor 


0.06 


4.74 


4.62 


4.46 








13141 


17848 


852 


02 49 01.5 


-62 48 23.5 


50.5 


A2V 


100 


1 




0.10 


5.25 


5.16 


4.97 








14551 


19545 


943 


03 07 50.8 


-27 49 52.1 


54.6 


A5V 


30 


3 


Tuc-Hor 


0.17 


6.18 


5.85 


5.77 








15648 


20677 


1002 


03 21 26.6 


+43 19 46.7 


46.2 


A3V 


125 


11 




0.05 


4.96 


4.86 


4.78 








16449 


21997 


1082 


03 31 53.6 


-25 36 50.9 


71.9 


A3IV/V 


30 


3 


Tuc-Hor 


0.12 


6.38 


6.12 


6.10 


yes 


1 




22192 


30422 


1525 


04 46 25.8 


-28 05 14.8 


56.2 


A3 


125 


11 




0.19 


6.18 


5.73 


5.72 


yes 


4 




22226 


30447 




04 46 49.5 


-26 18 08.8 


80.3 


F3V 


30 


2 


Columba 


0.39 


7.85 


6.95 


6.89 








23296 


32115 


1613 


05 00 39.8 


-02 03 57.7 


49.6 


A8 


125 


11 




0.29 


6.31 


5.62 


5.58 








26309 


37286 


1915 


05 36 10.3 


-28 42 28.9 


52.8 


A2 


30 


3 


Tuc-Hor 


0.15 


6.26 


5.94 


5.86 


yes 


3 




26624 


37594 


1940 


05 39 31.2 


-03 33 52.9 


42.6 


A8V 


125 


11 




0.29 


5.99 


5.36 


5.21 








32938 


50445 


2558 


06 51 42.4 


-36 13 49.0 


55.2 


A3V 


125 


11 




0.18 


5.94 


5.54 


5.51 








34782 


55568 


2720 


07 12 04.1 


-30 49 16.9 


47.6 


A8V 


125 


11 




0.27 


6.10 


5.46 


5.39 








35567 


57969 




07 20 23.0 


-56 17 40.7 


71.0 


A1V 


20 


9 




0.11 


6.57 


6.39 


6.28 








41152 


70313 


3277 


08 23 48.5 


+53 13 11.0 


50.4 


A3V 


200 


1 




0.13 


5.52 


5.29 


5.25 


yes 


1 




41307 


71155 


3314 


08 25 39.6 


-03 54 23.1 


37.5 


AOV 


400 


5 




-0.01 


3.91 


4.09 


4.08 


yes 


1 




42334 


73495 


3420 


08 37 52.2 


-26 15 18.0 


71.1 


A0V 


125 


11 




-0.03 


5.24 


5.35 


5.32 








44923 


78702 


3638 


09 09 04.2 


-18 19 42.8 


83.7 


A0/A1V 


200 


6 




0.00 


5.73 


5.64 


5.69 








53771 


95429 


4296 


11 00 08.3 


-51 49 04.1 


61.1 


A3 


125 


11 




0.19 


6.15 


5.79 


5.77 






yes/10" 


57013 


101615 


4502 


11 41 19.8 


-43 05 44.4 


65.5 


AOV 


125 


11 




0.04 


5.54 


5.51 


5.44 








57328 


102124 


4515 


11 45 17.0 


+08 15 29.2 


37.4 


A4V 


125 


11 




0.17 


4.84 


4.54 


4.41 








60595 


108107 


4722 


12 25 11.8 


-11 36 38.1 


70.5 


A1V 


125 


11 




0.04 


5.95 


5.91 


5.83 








61468 


109536 


4794 


12 35 45.5 


-41 01 19.0 


35.5 


A7V 


100 


7 




0.22 


5.12 


4.71 


4.57 








61960 


110411 


4828 


12 41 53.1 


+ 10 14 08.3 


36.3 


AOV 


100 


1 




0.08 


4.88 


4.76 


4.68 


yes 


1 




62983 


112131 


4901 


12 54 18.7 


-11 38 54.9 


68.5 


A2V 


125 


11 




0.07 


6.00 


5.85 


5.83 






yes/0.04" 


66634 


119024 


5142 


13 39 30.4 


+52 55 16.4 


53.6 


A3V 


125 


11 




0.11 


5.46 


5.20 


5.17 








69713 


125161 


5350 


14 16 09.9 


+51 22 02.0 


29.1 


A9V 


125 


11 




0.24 


4.75 


4.32 


4.29 






yes/32" 


69732 


125162 


5351 


14 16 23.0 


+46 05 17.9 


30.4 


AO 


200 


1 




0.09 


4.18 


4.03 


3.91 


yes 


1 




78078 


142703 


5930 


15 56 33.4 


-14 49 46.0 


51.1 


A2 


125 


11 




0.24 


6.11 


5.39 


5.34 








92024 


172555 


7012 


18 45 26.9 


-64 52 16.5 


28.5 


A7V 


12 


8 


P Pictoris 


0.20 


4.78 


4.25 


4.30 


yes 


1 


yes/72" 


95261 


181296 


7329 


19 22 51.2 


-54 25 26.2 


48.2 


AOV 


12 


8 


P Pictoris 


0.02 


5.02 


5.15 


5.01 








99273 


191089 




20 09 05.2 


-26 13 26.5 


52.2 


F5V 


12 


2 


P Pictoris 


0.48 


7.18 


6.12 


6.08 


yes 


1 




104365 


201184 


8087 


21 08 33.6 


-21 11 37.2 


55.1 


AOV 


125 


11 




0.00 


5.30 


5.33 


5.31 






yes/70" 


110935 


212728 


8547 


22 28 37.7 


-67 29 20.6 


43.1 


A4V 


125 


11 




0.21 


5.56 


5.14 


5.05 








114189 


218396 


8799 


23 07 28.7 


+21 08 03.3 


39.4 


A5V 


30 


3 


Columba 


0.26 


5.97 


5.28 


5.24 


yes 


1 




115738 


220825 


8911 


23 26 56.0 


+01 15 20.2 


47.1 


AO 


70 


3 


AB Dor 


0.04 


4.95 


4.95 


4.93 








116431 


221853 




23 35 36.2 


+08 22 57.4 


68.4 


FO 


100 


1 




0.41 


7.35 


6.44 


6.40 
















Additional targets included from 


previous surveys or discoveries. 
















27321" 


39060 


2020 


05 47 17.09 


-51 03 59.44 


19.4 


A6V 


12 


8 


P Pictoris 


0.17 


3.86 


3.54 


3.53 


yes 


1 




61498 6 


109573 


4796 


12 36 01.03 


-39 52 10.23 


72.8 


AO 


8 


10 


TWA 


0.00 


5.78 


5.79 


5.77 


yes 


1 


yes/8" 


98495* 


188228 


7590 


20 00 35.56 


-72 54 37.82 


32.2 


AOV 


50 


7 




-0.03 


3.95 


3.76 


3.80 









Notes. (£I) Source: Bonnefoy et al. (201 1). Data from ESO program 284.C-5057, PI A.-M. Lagrange. The data were drawn from the ESO archive and analyzed following the procedure described in 
Sect. 3. w Source: Chauvin et al. (2010). Detection limits provided by G. Chauvin (private communication). 



References. (1) Rheeetal. (2007); (2) Mooretal. (2011); (3) Zuckerman et al. (2011); (4) Su et al. (2006); (5) Janson et al. (2011); (6) Mooretal. (2006); (7) Chauvin et al. (2010); (8) 
Zuckerman et al. (2001); (9) Tetzlaff et al. (2011); (10) Staufferet al. (1995); (11) see Sect. 2. 
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to 5 pixel on the detector over the full observing sequence. For 
NaCo, a well-documented problem with the pupil-tracking mode 
(corrected in October 201 1) induced a slow drift of the star dur- 
ing some of the observations. This drift, which could reach sev- 
eral pixels per minute, is proportional to the star elevation. To 
avoid any important drift, the observing strategy included fre- 
quent recentering of the star during the observations and the use 
short DITs (during the follow-up phase of the survey). 

In each observing sequence, a frame selection was per- 
formed to reject the worst frames where a sudden increase of see- 
ing or an interruption of the AO loop occurred. The frame selec- 
tion was based on three parameters: the position of the PSF cen- 
ter in each frame compared to previous and subsequent frames, 
the overall flux in the frame and a measure of the fraction of flux 
contained in an annulus covering an area just outside of the sat- 
urated region. Frames for which a sudden change of the PSF po- 
sition or a drop in intensity are observed generally correspond to 
frames where the AO loop is open, and are thus of poor quality. 
These frames were inspected visually and in most cases removed 
from the observing sequence. After the frame selection process, 
the sequence of images usually contained frames of equivalent 
quality in terms of AO correction. 

For astrometric calibration of the NIRI data, we used archive 
data for the ©i Ori C field observed on November 30 2008 with 
the f/32 camera, K' filter and without the ALTAIR field lens. 
We used the coordinates reported by McCaughrean & Stauffer 
(1994) to determine the mean pixel scale and true North orien- 
tation of the NIRI detector. The mean pixel scale and true North 
orientation measured were respectively 21.38 + 0.14 mas/pixel 
and 0.27 + 0.05 deg. For NaCo observations, we used calibra- 
tion data collected within our observing programs and within 
the NaCo Large Program collaboration (NaCo-LP 1 ; Chauvin 
2010). The calibrators used were the ©i Ori C field and 
IDS 13022+0107, when the former was not observable. The cal- 
ibration was performed with the S13 camera and the H filter 
at several epochs between November 2009 and January 2012. 
The mean pixel scale measured was 13.19 + 0.07 mas/pixel, 
and the true North orientation varied from -0.17 ± 0.05 deg to 
-0.71 + 0.05 deg between the different epochs. The NaCo PSFs 
are usually not perfectly centro-symmetric, causing the centers 
determined on unsaturated and saturated PSFs to be slightly dif- 
ferent. The effect has been estimated to ~6 mas and constitutes 
the dominant astrometric error. It was estimated using a dedi- 
cated calibration where saturated and non-saturated PSFs were 
taken alternatively, and their center determined using Moffat or 
Gaussian profile fitting. 

Finally, the cleaned registered ADI images were processed 
with an implementation of the LOCI algorithm (Lafreniere et al. 
2007b). All the targets were analyzed using a similar set of pa- 
rameters: N s = 0.75 FWHM, N A = 300 for NaCo and N A = 500 
for NIRI. For each target, partial subtraction induced by LOCI 
was measured using fake planets introduced into the data at sev- 
eral position angles and angular separations, in a way similar to 
that described by Lafreniere et al. (2007b). Temporal smearing 
induced by the FoV rotation was also measured using fake plan- 
ets. A typical final reduced image is visible in Fig. 3 



1 ESO program 184.C-0567, PI J.-L. Beuzit, "Probing the 
Occurrence of Exoplanets and Brown Dwarfs at Wide Orbits". 




Fig. 3. Example of final image after data reduction and process- 
ing with LOCI for target HIP 10670. The central white region 
corresponds to the saturated area in the science images, which 
has been masked in the data analysis. The point source 4.7" West 
of the star is a background object 14.6 mag fainter than the cen- 
tral star. 



4. Results 

4.1. Identification of companion candidates 

Detection of the candidate companions (CCs) was performed by 
visual inspection of the final images after LOCI processing and 
of the signal-to-noise ratio (SNR) maps. The SNR maps were 
obtained by normalizing each pixel in the images to the noise 
measured in an annulus of width A/D at the same separation as 
the pixel from the star. The maps were used to identify CC with 
a detection level above 5<x. 

Relative astrometry for each CC was obtained using 2D 
Gaussian fitting, and the position was refined by maximizing the 
flux in a 1 .5 A/D aperture moved on a 2 x 2 pixel grid by steps 
of 0.1 pixel around the position found with the Gaussian fitting. 
The main astrometric error for the NaCo data is the 0.5 pixel un- 
certainty on the actual center of the PSF in saturated data (see 
Sect. 3). Other sources of error come from the shift induced by 
the use of the neutral density (0.1 pixel) or the 2D Gaussian 
fitting (usually 0.1 pixel). In NIRI data, we assume the domi- 
nant error term to be the residual distortion (6 mas) close to the 
star. At separations larger than 5", the distortion correction im- 
proves significantly the astrometric measurements for the CCs 
but not enough to reach a similar accuracy. For these candi- 
dates we have assumed a residual error of 1 pixel (21.4 mas). 
Although there are some uncertainties induced by the distortion 
in the NIRI field of view, the precision was sufficient to clearly 
identify background or comoving objects. An updated distortion 
map for NIRI will be presented in Galicheret al. (2012, in prepa- 
ration). 

Relative photometry was estimated by aperture photometry 
in 1.5A/D aperture centered on the best estimation of the posi- 
tion. Contrast difference with respect to the star was obtained by 
dividing the flux of each CC by the flux of the non-saturated PSF 
measured in an aperture of the same diameter, scaled by the ra- 
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Table 2. Summary of the observations. 



Name 


Date 


Instrument 


Filter 


Num. of exposures 


DIT 

(s) 


NDIT 


T 

A exp 

(min) 


FoV rotation 
(deg) 


Median seeing 
(as) 


Observed targets 


HIP 7345 


2007-09-12 


NIRI 


CH4s 


140 


30.0 




70 


38.9 


0.4 


HIP 10670 


2007-09-16 


NIRI 


CH4s 


140 


30.0 


[ 


70 


87.5 


0.5 




2008-10-12 


NIRI 


K' 


25 


30.0 


j 


13 


7.5 


0.5 


HIP 11360 


2007-12-30 


NIRI 


CH4s 


120 


30.0 


j 


60 


68.2 


0.5 




2008-10-12 


NIRI 


K' 


15 


30.0 




8 


6.2 


0.4 


HIP 12413 


2009-12-20 


NaCo 


K s 


125 


20.0 




42 


26.0 


1.1 




2011-11-07 


NaCo 


K s 


2200 


1.2 


j 


44 


33.3 


1.0 


HIP 13141 


2009-12-21 


NaCo 




125 


20.0 


j 


42 


15.7 


1.0 


HIP 14551 


2009-12-21 


NaCo 


K s 


125 


20.0 




42 


10.1 


1.2 


HIP 15648 


2008-10-17 


NIRI 


K' 


128 


30.0 


j 


64 


51.2 


0.7 


HIP 16449 


2007-09-22 


NIRI 


CH4s 


98 


30.0 


j 


49 


32.8 


0.9 




2008-10-12 


NIRI 


K' 


15 


30.0 


j 


8 


2.6 


0.4 


HIP 22192 


2009-12-20 


NaCo 


K s 


125 


20.0 




42 


22.4 


1.4 


HIP 22226 


2008-01-15 


NIRI 


CH4s 


120 


30.0 


j 


60 


33.4 


0.6 




2008-11-18 


NIRI 


K' 


26 


30.0 


j 


13 


4.8 


0.8 


HIP 23296 


2010-01-03 


NIRI 


K' 


252 


15.0 


j 


63 


56.7 


0.5 


HIP 26309 


2009-12-21 


NaCo 


K s 


125 


20.0 




42 


102.1 


0.8 


HIP 26624 


2008-11-14 


NIRI 


K' 


160 


30.0 


j 


80 


55.8 


0.4 


HIP 32938 


2009-12-20 


NaCo 


K s 


125 


20.0 


j 


42 


45.9 


0.8 


HIP 34782 


2009-12-21 


NaCo 




125 


20.0 


j 


42 


77.4 


0.8 




2012-01-11 


NaCo 


K s 


1000 


2.8 




47 


77.9 


0.9 


HIP 35567 


2009-12-20 


NaCo 


K s 


100 


20.0 


j 


33 


17.7 


1.4 




2011-12-18 


NaCo 


K s 


360 


7.5 


j 


45 


20.9 


0.8 


HIP 41 152 


2010-01-03 


NIRI 


K' 


232 


15.0 


j 


58 


38.2 


0.6 


HIP 41307 


2009-12-21 


NaCo 


K s 


100 


20.0 


j 


33 


17.9 


1.1 


HIP 42334 


2009-12-20 


NaCo 


K s 


125 


20.0 


j 


42 


19.6 


1.4 




2011-04-27 


NaCo 




530 


2.5 


j 


22 


5.5 


1.1 


HIP 44923 


2008-03-23 


NIRI 


CH4s 


120 


30.0 




60 


40.8 


0.5 


HIP 53771 


2010-03-06 


NaCo 


K s 


25 


20.0 




42 


19.6 


1.0 




2012-01-12 


NaCo 


K s 


880 


4.0 




59 


29.1 


0.8 


HIP 57013 


2010-03-07 


NaCo 


K s 


25 


20.0 




42 


29.2 


0.6 


HIP 57328 


2008-03-24 


NIRI 


CH4s 


130 


30.0 


j 


65 


97.8 


0.4 


HIP 60595 


2010-03-07 


NaCo 


K s 


25 


20.0 




42 


30.1 


0.5 


HIP 61468 


2009-06-25 


NaCo 


K s 


125 


20.0 




42 


33.4 


1.0 


HIP 61960 


2008-03-21 


NIRI 


CH4s 


122 


30.0 




61 


99.8 


0.4 


HIP 62983 


2010-03-06 


NaCo 


K s 


25 


20.0 




42 


37.8 


0.9 


HIP 66634 


2009-02-03 


NIRI 


K' 


42 


30.0 


j 


21 


10.7 


0.7 


HIP 69713 


2008-03-21 


NIRI 


CH4s 


140 


30.0 


j 


70 


51.2 


0.5 


HIP 69732 


2009-02-02 


NIRI 


K' 


88 


30.0 


j 


44 


33.9 


0.5 


HIP 78078 


2009-06-24 


NaCo 


K, 


100 


20.0 


j 


33 


40.7 


1.4 


HIP 92024 


2009-06-25 


NaCo 


K s 


230 


20.0 


j 


77 


32.0 


1.1 


HR 7329 


2008-08-07 


NaCo 


H 


27 


25.0 




23 


11.0 


0.5 


HIP 99273 


2007-08-10 


NIRI 


CH4s 


140 


30.0 


"[ 


70 


34.4 


0.5 




2008-09-04 


NIRI 


CH4s 


130 


30.0 


1 


65 


32.2 


0.5 


HIP 104365 


2007-09-18 


NIRI 


CH4s 


140 


30.0 




70 


45.4 


0.5 




2008-09-09 


NIRI 


CH4s 


20 


30.0 




10 


5.2 


0.8 


HIP 110935 


2009-06-24 


NaCo 




125 


20.0 




42 


14.4 


1.6 


HR 8799 


2007-10-17 


NIRI 


CH4s 


120 


30.0 




60 


172.3 


0.4 




2008-09-01 


NIRI 


CH4s 


70 


30.0 




35 


152.4 


0.6 




2008-09-05 


NIRI 


K' 


60 


30.0 




30 


139.1 


0.3 




2008-10-10 


NIRI 


CH4s 


120 


30.0 




60 


171.6 


0.6 




2008-10-14 


NIRI 


K' 


60 


12.0 




12 


131.2 


0.4 


HIP 115738 


2008-10-15 


NIRI 


CH4s 


120 


30.0 




60 


69.4 


0.4 


HIP 116431 


2007-09-12 


NIRI 


CH4s 


140 


30.0 




70 


72.3 


0.4 




2008-09-24 


NIRI 


CH4s 


60 


30.0 




30 


54.9 


0.7 


Targets from previous imaging surveys 


P Pictoris 


2010-04-10 


NaCo 




160 


0.15 


150 


60 


21.0 


0.8 


HIP 61498 


2003-06-06 


NaCo 


H 


12 


3.0 


25 


15 




0.7 


HIP 98495 


2004-04-27 


NaCo 


H 




3.0 


20 






0.6 
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tio of exposure time and the transmission of the neutral density 
(NaCo data) or the correction due to the fact that unsaturated ob- 
servations were not performed in the same filters as saturated ob- 
servations (NIRI data). This correction factor for NIRI was cal- 
culated using slightly saturated images acquired regularly during 
the ADI sequence with the same filter as the deeply saturated 
science images. The non-saturated PSF was normalized to these 
slightly saturated images using the integrated flux in an annulus 
outside of the saturated part of the images (usually between 10 
and 20 pixels of radius). A calibration performed on data satu- 
rated and unsaturated with the same filters indicated an error on 
the flux normalization of ~4%. The flux loss induced by LOCI 
partial subtraction and by the temporal smearing of the CC PSF 
by field rotation was also included as a correction term. Finally, 
the photometric error was estimated by including a measure of 
the local SNR at the separation of each CC and the error on the 
calibration of the flux normalization between the saturated and 
unsaturated PSF. 

For newly detected stellar companions (see Sect. 4.3.1), 
where a more robust estimation of the relative astrometry and 
photometry is useful for the determination of their physical pa- 
rameters, a different approach was taken: instead of estimating 
the position and signal level of the CC in the final reduced im- 
age, we implemented a routine to subtract a "negative CC" in 
each raw frames before performing the ADI optimization in a 
way similar to that described in Marois et al. (2010a). The rou- 
tine was run multiple times for various insertion position of the 
negative CC in order to minimize the residual noise in the final 
image at the expected position of the CC and in the surrounding 
area. When a position minimizing the residuals was found, the 
routine was used again for different injection flux for the CC, 
again trying to minimize the residual noise in the final image. 
Errors were estimated using variations of the best injection point 
and flux obtained by measuring the residuals in areas of different 
sizes around the position of the CC. 

A total of ~50 CCs were identified above the 5cr level around 
22 of the observed targets. Targets with CCs at a projected phys- 
ical separation less than or equal to 320 AU were followed- 
up by subsequent observations to verify if they are co-moving 
with the target star. This cut-off was chosen according to the 
widest likely known companion at the time of the follow-up, 
1RXS J160929. 1-210524 b, which is at a projected separation 
of ~330 AU (Lafreniere et al. 2008). There is only one target 
for which no second epoch was obtained although there was one 
CC within 320 AU: HIP 41307 had a CC at a projected sep- 
aration of 217 AU but it was confirmed as a background star 
by Janson et al. (201 1). For targets with follow-up observations, 
the position of the candidates was checked for consistency with 
either a background or a co-moving object based on the target 
parallax and proper motion, also taking into account uncertain- 
ties on these particular values. Similarly to Chauvin et al. (2010) 
and Ehrenreich et al. (2010), we performed ax 2 probability test 
on the difference in position Aa and AS with respect to the star 
at two epochs. These probability tests along with the proper mo- 
tion plots in RA/DEC or in separation/position angle for each 
CC allowed attributing a status of background, comoving or am- 
biguous to each of them. The properties of all detected CCs are 
summarised in Appendix A. 

4.2. Detection limits 

The ultimate contrast reached by our observations is quantified 
using 5cr detection limits obtained by measuring at each angu- 
lar separation the residual noise in an annulus of width A/D, 
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Fig. 4. Examples of /f-band 5cr detection limits representing the 
best (HIP 69732) and worst (HIP 15648) limits, as well as an in- 
termediate case (HIP 26624). The right axis shows the detection 
limit in masses of Jupiter for the median target of our survey, 
a 100 Myr-old A3V star at 50 pc, converted using the COND 
evolutionary models (Baraffe et al. 2003). 



normalized by the flux of the unsaturated PSF. As mentioned 
previously, two correction factors were applied to the detection 
limits. The first one accounts for the flux loss of point sources 
induced by the LOCI processing. It is calculated by introduc- 
ing fake companions of known flux at various separation and 
position angle into the raw data, fully processing the data with 
LOCI, and measuring the flux of the companions in the final im- 
age (see Lafreniere et al. 2007b). This process was repeated 10 
times with the companions at various azimuthal positions to ob- 
tain the average attenuation as a function of angular separation 
from the star. The second correction factor takes into account 
the smearing of the PSF of an object in the field induced by the 
field rotation during the exposures. The dilution factor was cal- 
culated by simulating the smearing of point sources at various 
separations for every frame. Finally, the detection limits were 
corrected by dividing them at all separations by the appropriate 
attenuation for the flux loss and dilution for the smearing. 

The best magnitude difference detection limits for angular 
separations ranging from 0.5" to 10.0" are listed in Table 3. 
Example limits representative of good, median and poor contrast 
performance in /f-band are plotted in Fig. 4. The observations 
are background-limited outside ~5" for all targets. The back- 
ground limits range between 19 and 23 mag for NIRI in //-band, 
between 17.5 and 22.5 mag for NIRI in Zf-band, and between 
16.5 and 20.5 mag for NaCo in Zf-band. These values are within 
the typical ranges for observations in good conditions with NIRI 
and NaCo in 0.5 to 1 hour of integration time. 

A common concern for the estimation of detection lim- 
its in high-contrast data is that of the detection threshold. It 
is known that speckle noise does not follow Gaussian statis- 
tics (Goodman 1968; Fitzgerald & Graham 2006). However, 
Marois et al. (2008a) have performed a detailed analysis of NIRI 
data processed with ADI (Marois et al. 2006) and demonstrated 
that the residual noise follows quasi-Gaussian statistics when a 
sufficient number of images (i.e. of FoV rotation) are combined. 
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Table 3. Best 5cr detection limit for each target of the sample. 



Name 


0.50" 


0.60" 


0.75" 


1.00" 


2.00" 


3.00" 


4.00" 


5.00" 


7.50" 


10.00" 


Observed targets 


HIP 7345 


10.6 


11.5 


12.3 


13.3 


15.4 


16.3 


16.8 


16.9 


17.1 


17.1 


HIP 10670 








13.5 


15.8 


16.7 


17.2 


17.6 


18.0 


18.1 


HIP 11360 


10.6 


11.6 


12.2 


13.1 


15.1 


16.0 


16.4 


16.6 


16.7 


16.7 


HIP 12413 


8.9 


10.0 


11.4 


12.6 


14.2 


14.5 


14.6 


14.6 


13.6 




HIP 13141 


9.2 


10.0 


11.3 


12.7 


14.6 


15.1 


15.3 


15.4 


14.8 


14.4 


HIP 14551 


8.5 


9.3 


10.8 


12.3 


14.3 


14.8 


14.9 


15.0 


14.4 


14.2 


HIP 15648 




9.1 


10.1 


10.9 


12.3 


12.7 


12.8 


12.8 


12.9 


12.8 


HIP 16449 


8.0 


8.9 


9.6 


10.6 


12.5 


12.9 


13.0 


13.1 


13.2 


13.2 


HIP 22192 


9.5 


10.2 


11.5 


12.6 


14.2 


14.4 


14.5 


14.5 


13.8 


13.6 


HIP 22226 


9.8 


10.7 


11.4 


12.3 


14.3 


15.0 


15.2 


15.4 


15.4 


15.4 


HIP 23296 




11.0 


11.9 


13.2 


15.4 


16.4 


16.7 


16.9 


16.9 


16.9 


HIP 26309 


10.0 


10.8 


11.9 


12.6 


14.1 


14.4 


14.4 


14.3 


13.2 


12.5 


HIP 26624 








11.8 


14.2 


15.3 


15.7 


15.8 


15.9 


15.9 


HIP 32938 


10.7 


11.6 


12.5 


13.3 


14.4 


14.5 


14.6 


14.6 


13.7 


12.8 


HIP 34782 


10.9 


11.8 


12.7 


13.8 


15.0 


15.3 


15.3 


15.3 


14.3 


13.4 


HIP 35567 


9.5 


10.9 


12.2 


13.3 


14.2 


14.4 


14.3 


14.4 


13.5 




HIP 41 152 






11.2 


12.8 


15.0 


16.2 


16.6 


16.8 


17.0 


17.0 


HIP 41307 




10.5 


11.7 


13.1 


15.4 


16.1 


16.4 


16.5 


16.1 


15.7 


HIP 42334 


9.7 


10.7 


11.5 


12.7 


14.4 


14.7 


14.9 


14.9 


14.2 


13.9 


HIP 44923 




10.4 


11.1 


11.8 


14.2 


15.4 


16.0 


16.2 


16.5 


16.5 


HIP 53771 


10.0 


11.2 


12.5 


13.6 


14.7 


14.8 


14.8 


14.8 


14.0 




HIP 57013 


8.7 


9.7 


10.8 


11.8 


13.8 


14.4 


14.6 


14.7 


14.4 


14.0 


HIP 57328 








13.5 


15.5 


16.7 


17.4 


17.7 


18.1 


18.2 


HIP 60595 


8.9 


9.6 


10.7 


11.8 


13.8 


14.3 


14.6 


14.7 


14.0 


13.5 


HIP 61468 




10.1 


11.2 


12.4 


14.5 


15.1 


15.2 


15.3 


14.4 


15.7 


HIP 61960 








12.7 


14.8 


16.0 


16.7 


17.0 


17.4 


17.4 


HIP 62983 


8.6 


9.4 


10.5 


11.6 


13.5 


14.0 


14.1 


14.2 


13.4 


13.3 


HIP 66634 




9.1 


10.3 


11.3 


13.6 


14.5 


14.7 


14.8 


14.8 


14.8 


HIP 69713 








12.4 


14.7 


15.8 


16.6 


16.9 


17.4 


17.5 


HIP 69732 










16.3 


17.4 


17.9 


18.2 


18.4 


18.5 


HIP 78078 


7.9 


9.0 


10.2 


11.2 


12.8 


13.2 


13.3 


13.3 


12.3 




HIP 92024 






10.8 


11.9 


13.9 


14.5 


14.6 


14.7 


9.7 




HR 7329 




11.2 


12.1 


13.1 


15.1 


16.0 


16.5 


16.7 






HIP 99273 


10.2 


11.3 


12.2 


13.4 


15.4 


16.1 


16.4 


16.6 


16.7 


16.7 


HIP 104365 




11.5 


12.4 


13.6 


15.5 


16.3 


16.9 


17.1 


17.3 


17.4 


HIP 110935 


8.5 


9.0 


9.9 


11.2 


13.1 


13.5 


13.6 


13.7 


13.2 




HR 8799 




12.0 


12.7 


13.6 


15.6 


16.4 


16.8 


17.0 


17.1 


17.0 


HIP 115738 




12.0 


12.9 


13.5 


15.6 


16.6 


17.1 


17.3 


17.5 


17.5 


HIP 116431 


10.7 


12.0 


12.9 


13.9 


15.8 


16.3 


16.5 


16.6 


16.7 


16.6 


Targets from previous imaging surveys 


P Pictoris 


10.3 


11.2 


12.0 


13.5 


15.2 


15.3 










HIP 61498 


9.5 


10.1 


10.4 


11.1 


13.4 


13.6 


14.0 


14.1 






HIP 98495 


10.6 


10.6 


10.7 


13.2 


14.7 


14.6 


14.6 


14.9 







Notes. Values are given in magnitude difference in the filter of the observations: ACH4s or AK' for NIRI observations, AA" S for NaCo observations. 



For most targets-, we obtained a FoV rotation of at least 20 deg 
(Table 2) over several tens of images, which should be sufficient 
to reach quasi-Gaussian noise statistics after LOCI processing. 
For comparison, previous surveys using some form of ADI ob- 
servations have adopted 5<x or 6cr thresholds to derive their de- 
tection limits, either assuming Gaussian statistics for the residual 
noise (Biller et al. 2007; Lafreniere et al. 2007a; Leconte et al. 
2010; Janson et al. 201 1) or using a higher threshold to account 
for the non-Gaussianity of the noise (Kasper et al. 2007). 

The 5<t detection limits were finally converted to absolute 
magnitudes using the Hipparcos parallax and 2MASS H or K s 
apparent magnitudes. Since observations were obtained in a fil- 
ter set different from the 2MASS filters (in particular the CH4s 

2 The second epochs with NIRI have usually been obtained with 
much shorter integration times than the first epochs, providing a much 
smaller FoV rotation. They were only intended to confirm the CCs de- 
tected in the first epoch images. 



NIRI filter), we checked that for the different spectral types of 
our targets (from A0V to F5 V) there was no significant change in 
apparent magnitude between the 2MASS filters and the NaCo or 
NIRI filters. The correction terms were found to be completely 
negligible (< 0.03 mag) compared to the other systematic uncer- 
tainties in the photometry. 

4.3. Companion detections 

Within the 42-star sample, 5 stars have newly identified or previ- 
ously reported physically associated companions, and we briefly 
discuss each case of a stellar or substellar companion in this 
section. Of the companion candidates with follow-up newly re- 
ported in this study, three are confirmed as co-moving, and all 
but one of the remaining candidates have a <0.01% probabil- 
ity to be comoving. The possible residuals on astrometry are 
well below the threshold that would change the assessment from 
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Fig. 5. Offset position (black circles) with uncertainties measured relative to the primary star for the co-moving stellar companions 
to HIP 42334 (left) and HIP 104365 (right). The variation expected for a background stationary object is also plotted (grey line) 
with the expected positions of background objects (grey square). 



background to co-moving. In one case, one of candidate com- 
panions to HIP 10670, it cannot be excluded that this object has 
a proper motion of its own, making it inconsistent with both a 
stationary background object and a co-moving object. Because 
it shows a motion much larger than a pixel (13.19 mas on NaCo, 
21.38 mas on NIRI) between the two epochs, suggesting a sig- 
nificant motion not compatible with a co-moving object, it is 
treated as a background object in the statistics. 



4.3.1. Stellar companions 

We have identified three co-moving stellar companions around 
two stars in our data. Both stars were already known to be in 
wide multiple systems, but the companions we report were pre- 
viously unknown due to their large contrast ratio with their par- 
ent star. Proper motion plots are given in Fig. 5, and relative 
positions and fluxes are summarized in Table 4. 

HIP 42334 - The star is a 125 Myr old A0V star at 71.07 pc 
in a double system, with the secondary component at a sepa- 
ration of 20" and a magnitude of 13 in the optical (Mason et al. 
2001). A point source 8.2+0.1 mag fainter in K s than the primary 
was detected at a separation of 3.176" + 0.007" and position an- 
gle (PA) of 1 1.52° + 0.12°. The star was re-observed because of 
the presence of other fainter point sources in the field, but only 
the brightest candidate proved to be co-moving (Fig. 5, left). 
Comparison to evolutionary models at 125 Myr (Chabrieret al. 
2000) place this object at the bottom of the main sequence, con- 
sistent with an M6-M8 dwarf. 

HIP 104365 - The star is 125 Myr old A0V star at 55.12 pc 
in a quadruple system (Mason et al. 2001). A very tight pair of 
companions was detected at a separation of ~9.57" and a PA 
of ~14.8°. After one year they both clearly showed a common 
proper motion with HIP 104365 (Fig. 5, right). Given their pro- 
jected separation of 9.8 AU and their common proper motion, it 
is extremely likely that they are gravitationally bound together. 



Table 5. Properties of the substellar companions from the sam- 
ple. 



Name 


Orbital sep. 


Mass 


References 




(AU) 


(M Jup ) 




HR 8799 b 


-68° 


c+2.0 
-L5 


1,2,3,4 


HR 8799 c 


-38" 


7 +3 
'-2 


1,2,3,4 


HR 8799 d 


-24" 


1 I 


1,2,3,4 


HR 8799 e 


-15" 


ii 


1,2,3,4 


ft Pictoris b 


8-9 


7-12 


5,6,7,8,9 


HR 7329 B 


-200" 


20-50 


10 



Notes. (n) Projected separation. 

References. (1) Marois et al. (2008a); (2) Marois et al. (2010b); (3) 
Currieetal. (2011); (4) Soummeretal. (2011); (5) Lagrange et al. 
(2009); (6) Lagrange et al. (2010); (7) Bonnefoy et al. (2011); (8) 
Lagrange et al. (2012a); (9) Chauvin et al. (2012); (10) Neuhauser et al. 
(2011). 



Their contrast difference with the star in the CH4s filter makes 
them compatible with M5-M7 dwarfs. 

4.3.2. Substellar companions 

Three targets in the sample - HR 8799, f3 Pic, and HR 7329 
- have previously reported substellar companions (Marois et al. 
2008b, 2010b; Lagrange et al. 2009; Lowrance et al. 2000; 
Neuhauser et al. 2011) including four planets and one brown 
dwarf. A summary of these objects and their properties is given 
in Table 5. We note that the observations of these three targets 
were not unusually sensitive and biased, but have detection lim- 
its similar to many other targets in the sample. 

HR 8799 - This is the first imaged multiple planet sys- 
tem, with four known planetary-mass companions orbiting the 
star at 68, 38, 24 and 14.5 AU (Marois et al. 2008b, 2010b). 
It is an A5V star located at 39.4 + 1.0 pc from the Sun, clas- 
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Table 4. Properties of the detected stellar companions. 



Primary name 


Separation 


PA 


ACH4s 




Mass" 




(as) 


(deg) 


(mag) 


(mag) 


(M G ) 


HIP 42334 


3.176 + 0.007 


11.52 ±0.12 




8.2 ±0.1 


0.07-0.11 


HIP 104365 


9.545 ± 0.006 


15.27 ± 0.04 


7.4 ±0.1 




0.09-0.13 




9.595 ± 0.006 


14.27 ± 0.04 


7.4 ±0.1 




0.09-0.14 



Notes. (n> The mass range was estimated using the DUSTY evolutionary models (Chabrier et al. 2000) taking into account the uncertainties on the 
distance, on the magnitude of the primary, on the contrast of the companion, and on the age of the system. 



sifted as y Doradus and A Bootis. It shows a far-infrared ex- 
cess emission (Vega-like) consistent with the presence of a warm 
dust disk within 6-15 AU, a massive cold dust disk within 90- 
300 AU and a dust halo extending up to 1000 AU (Rhee et al. 
2007; Su et al. 2009). Recent estimations of the system age 
converge towards 30 Myr (Doyon et al. 2010; Zuckerman et al. 
201 1; Currie et al. 201 1), with HR 8799 being a likely member 
of the Columba Association. This estimation leads to a mass of 
5-10 M Jup for HR 8799 cde, and 3.5-7 M Jup for HR 8799 b. 
This system has been extensively studied in the literature since 
its discovery and we do not intend to make an exhaustive re- 
view of the numerous results. For detailed studies, see for ex- 
ample Bowler et al. (2010b), Currie et al. (201 1), Barman et al. 
(201 la), Galicher et al. (201 1) and Soummer et al. (201 1). 

j3 Pictoris - This Myr-old (Zuckerman et al. 2001) 

A5V star at a distance of 19.4 + 0.2 pc (Crifo et al. 1997) 
has been well-known for two decades since the imaging of a 
debris disk (Smith & Terrile 1984) with several asymmetries 
and a warp at -80 AU (Kalas & Jewitt 1995; Mouillet et al. 
1997), which could be explained by the presence of a gi- 
ant planet within 50 AU of the star. Lagrange et al. (2009, 
2010) confirmed the presence of a planetary-mass compan- 
ion with an orbital separation compatible with formation by 
core-accretion. Using the COND evolutionary models and K- 
band observations, Bonnefoy et al. (2011) infered a mass range 
of 7-10 Mj up for the planet. Follow-up of the system with 
imaging (Lagrange et al. 2012a; Chauvin et al. 2012) and RV 
measurements (Lagrange et al. 2012b) recently allowed an ac- 
curate determination of the orbital parameters and the mass 
of the planet. Using Markov Chain Monte-Carlo methods, 
Chauvin et al. (2012) measure a probable range of a = 8-15 AU 
for the semimajor axis, with a low eccentricity e < 0.16 and a 
high inclination i = 88.5 + 1.7 deg. Using a dataset where the 
planet, the main disk and the warp are detected, (Lagrange et al. 
2012a) also confirmed that the planet is located above the mid- 
plane of the main disk, making it aligned with the warp. These 
observations strongly support that the presence of the planet 
is responsible for the warped morphology of the disk. Finally, 
using 8 years of RV measurements on the star, Lagrange et al. 
(2012b) determine an upper limit of 10-12 Mj up for the mass of 
the planet, assuming a semimajor axis of 8-9 AU. 

HR 7329 - The substellar companion to the f3 Pic mov- 
ing group member HR 7329 (Zuckerman et al. 2001) was dis- 
covered with coronagraphic imaging with the Hubble Space 
Telescope (Lowrance et al. 2000). It has a magnitude difference 
of AF160W = 6.9 + 0.1 and separation of 4.170" + 0.005" 
(200 AU at 47.7 pc). It was further characterized with spectro- 
scopic observations in the //-band (Guenther et al. 2001). The 
most comprehensive compilation of astrometric and photometric 
measurements of this system and a high significance confirma- 
tion that the companion is physically associated with the primary 



is given in Neuhauser et al. (201 1). Like the HR 8799 and (3 Pic 
systems, HR 7329 has a debris disk (Smith et al. 2009). 

5. Statistical analysis 

5.1. Context 

Using AO instruments on large ground-based telescopes, a 
number of deep imaging surveys have been carried out to 
search for massive planets around young nearby solar-type stars 
(Chauvin et al. 2003; Neuhauser et al. 2003; Masciadri et al. 
2005; Billeretal. 2007; Kasperetal. 2007; Lafreniere et al. 
2007a; Chauvin et al. 2010; Delorme et al. 2012), and more re- 
cently around much more massive stars (Jansonetal. 2011). 
Despite a large improvement of sensitivity to massive planets in 
the recent years, due largely to the development of advanced data 
analysis methods to suppress the speckle noise, none of these 
surveys have reported the detection of planetary objects poten- 
tially formed by core-accretion, i.e. within the circumstellar disk 
of the primary. 

Taking advantage of their null results, the most recent large- 
scale surveys have placed constraints on the population of long- 
period massive planets. Using a combination of Monte-Carlo 
simulations and a Bayesian approach, it is possible to estimate 
the upper limit of the fraction of stars having substellar compan- 
ions (/ max ) within a certain range of orbital separation and mass. 
Nielsen & Close (2010) performed a detailed statistical analysis 
using the data for 118 stars from the surveys of Masciadri et al. 
(2005), Biller et al. (2007) and Lafreniere et al. (2007a), which 
allowed them to constrain the presence of planets down to 4 Mj up 
around F-M stars. Assuming a flat distribution for the mass and 
semimajor axis of planets, they rule out at 95% confidence the 
presence of planets more massive than 4 Mj up around more than 
20% of stars in the range 80-280 AU using the cold-start evo- 
lutionary models (see Sect. 5.3 for details on the evolutionary 
models) and in the range 20-500 AU using the hot-start mod- 
els, which is in agreement with the values previously reported 
by Lafreniere et al. (2007a). From the observations of 22 very 
young G-M targets in the Tucana and /3 Pic moving groups, 
Kasper et al. (2007) place an upper limit of 5% for planets above 
2-3 Mj up at separations larger than 30 AU. Using a sample of 
88 young nearby stars, Chauvin et al. (2010) report with 95% 
confidence that less than 10% of stars can have planetary-mass 
companions beyond 40 AU, assuming mass and semimajor axis 
distributions similar to that reported by Cumming et al. (2008) 
from RV surveys. Finally, Janson et al. (201 1) compared predic- 
tions of disk instability models to the observation of 15 B2-A0 
stars, and they showed that less than 30% of massive stars can 
form and retain companions below 100 Mj up within 300 AU at 
99% confidence. 

Previous surveys have also used their non-detections to in- 
vestigate with what confidence planet distributions obtained 
from RV surveys can be extrapolated for planets at much larger 
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orbital radii. This is in particular the case for the mass and period 
and semimajor axis distributions, which appear to be well fit- 
ted by simple power laws (Cumming et al. 2008) up to a certain 
semimajor axis cutoff. For the power law and planet frequency 
values reported by (Cumming et al. 2008), Nielsen & Close 
(2010) place the semimajor axis cutoff within 180 AU at 95% 
confidence with the cold-start models, and within 65 AU with 
the hot-start models. Similarly, Kasper et al. (2007) rule out 
the possibility of a positive power law index for a semima- 
jor axis distribution with a cutoff beyond 30 AU. With their 
data, Chauvin et al. (2010) explored a wide range of parameters 
for the power law indices and concluded that the non-detection 
probability is very sensitive to the semimajor power law index, 
which means that strong constraints on this parameter can be 
inferred from observations. For the combinations of parameters 
reported by Cumming et al. (2008), Chauvin et al. (2010) report 
a value of the semimajor axis cutoff within -100 AU with 95% 
confidence. 

Although the reported values are not identical due to the dif- 
ferent samples and sensitivities, all previous surveys converge 
toward a similar conclusion, which is that the frequency of plan- 
ets above 2-4 Mj up around F-M stars is below ~ 10-20% beyond 
~50 AU at high-confidence level. They also show that the distri- 
bution of semimajor axis is consistent with a power law of neg- 
ative index, i.e. with a decrease of the number of planets with 
orbital separation. 

The present work is similar in its aims to previous surveys, 
but with two major exceptions: it is focused on A-stars, and 
for the first time it includes the detections of the planetary sys- 
tem around HR 8799 and the planet around /3 Pic. A-stars are 
particularly interesting because RV surveys show a strong cor- 
relation between planet mass and stellar mass (Johnson et al. 
2007, 2010a). Although this correlation may not hold for plan- 
ets at large radii, it is tempting to test this hypothesis, espe- 
cially given the planets recently imaged around young A-type 
stars (Marois et al. 2008b, 2010b; Lagrange et al. 2009, 2010). 
It is therefore essential that our statistical analysis take into ac- 
count the confirmed detections for the estimation of the planet 
frequency around A-stars. 

In the following sections, we present our statistical formal- 
ism and several applications of our data to place constraints on 
the frequency of giant exoplanetary systems at large orbital radii 
around A-stars, and their possible mass and semimajor axis dis- 
tributions. We perform the analysis for a total of 42 stars: the 
39 that we have observed with NaCo and NIRI, and 3 young 
A-stars included from the survey of young, nearby austral stars 
by Chauvin et al. (2010) - /? Pic (using the data presented in 
Bonnefoy et al. 201 1), HIP 61498 and HIP 98495. 

5.2. Statistical formalism 

Our formalism for the statistical analysis is based on previous 
works by Carson et al. (2006) and Lafreniere et al. (2007a). We 
consider the observation of N stars enumerated by j = 1 . . .N. 
We note / the fraction of stars that have at least one compan- 
ion of mass and semimajor axis in the interval [m m i a ,m max ] n 
[flmin, a max ], and pj the probability that such a companion would 
be detected from our observations. The analysis does not differ- 
entiate between a single planet system and a multiple planet sys- 
tem, i.e. / represents the frequency of giant planetary systems at 
large orbital radii. With these notations, the probability of detect- 
ing such a companion around star j is (fpj) and the probability 
of not detecting it is (1 - fpj). Denoting {dj} the detections made 
by the observations, such that dj equals 1 if at least one compan- 



ion is detected around star / and otherwise, the likelihood of 
the data given / is 



Lad J w = Y](i-fpj) 1 - dj -(fPj) d 



(1) 



7=1 



The determination of the probability that the fraction of stars 
having at least one companion is / is obtained from Bayes' the- 



....... L({dj}\f)-p(f) 

Piflidj}) = -j , 

X L{{dj}\f) ■ p(f)df 



(2) 



where p(f) is the a priori probability density of /, or prior distri- 
bution, and p{f\{dj}) is the probability density of / given the ob- 
servations {dj}, or posterior distribution. Since we have no a pri- 
ori knowledge of the wide-orbit massive planets frequency, we 
adopt a "maximum ignorance" prior, p(f) = 1 . This approach is 
similar to previous studies for deep imaging surveys, allowing a 
direct comparison of the results. 

Given a confidence level (CL) a, we can use the posterior 
distribution p(f\{dj}) to determine a confidence interval (CI) for 
/ using 



P (f\{dj))df, 



(3) 



where / m ; n and / max represent the bounding values of /, i.e. the 
minimal and maximal fraction of stars with at least one planetary 
companion. Previous studies have always assumed a value of 
/min = since no detection was included. In that case, Eq. 3 be- 
comes an implicit equation on / max , which can be solved numer- 
ically to estimate / max . In our case, we want to use our detections 
to infer a value for both f^ m and / max . Following Lafreniere et al. 
(2007a), an equal-tail CI [/min,/m ax ] is found by solving: 



l-a _ f 1 

~ L 

1 - a C 1 ' 

= Jo 



P(f\{dj})df 



P(f\{dj})df 



(4) 
(5) 



Again, when a value of a is fixed, Eq. 4 and 5 become implicit 
equations on / m ; n and / max that can be solved numerically. 

5.3. Monte-Carlo simulations 

The critical step of the statistical analysis is the determination 
of pj, the completeness, for all observed targets. This complete- 
ness represents the fraction of companions that would have been 
detected in the interval [m m j n , f« max ] n [flmin, a max ] with the obser- 
vations. It is directly related to the sensitivity of the observations 
to the physical parameters of the planets, which in turn is a func- 
tion of the detection limits, age and distance of each star. 

A classical approach to estimate pj is to use Monte-Carlo 
(MC) simulations to generate large populations of planets with 
random physical and orbital parameters, and check their de- 
tectability given the sensitivity limits of the survey. For the MC 
simulations, we used the Multi-purpose Exoplanet Simulation 
System (MESS, Bonavita et al. 2012), a tool specifically de- 
signed for the statistical analysis of exoplanet surveys, which has 
already been used for the statistical analysis of several surveys 
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HR 8799 - 30 Myr - 39.4 pc HR 8799 - 30 Myr - 39.4 pc 




Mass (M Jup ) Mass (M jLlp ) 

Fig. 6. The detection limit (dashed line) for wide orbit planets from the observations of the target HR 8799 in the NIRI CH4s filter 
(left) and K' filter (right) compared to the DUSTY (red line) and COND (black line) evolutionary models (Chabrier et al. 2000; 
Baraffe et al. 2003). The DUSTY model mass-magnitude grid does not extend to magnitudes as faint as the detection limits. The 
current estimates of the HR 8799 planet masses (Marois et al. 2010b) are also plotted, and the empirical masses are consistent with 
the COND magnitudes within the error bars. 



(Chauvin et al. 2010; Janson et al. 2011; Delorme et al. 2012). 
The mass and semimajor axis of the planets are generated ei- 
ther from a linear grid or from specific distributions like power 
laws extrapolated from RV surveys of close-in planets (e.g. 
Fischer & Valenti 2005; Cumming et al. 2008; Johnson et al. 
2010a). Then for each combination of mass/semimajor axis, a 
large number of planets are generated by sampling randomly the 
other orbital parameters: the inclination has a uniform distribu- 
tion in sin i, the longitude of the ascending node and the argu- 
ment of periastron have a uniform distribution between and 
2n, and the time of passage at periastron has a uniform distribu- 
tion between and P, with P the orbital period. The eccentric- 
ity distribution is uniform between < e < 0.8 as suggested 
by surveys of long period exoplanets (Cumming et al. 2008). 
However, given the large number of planets generated at each 
mass/semimajor axis point, the effect of the eccentricity distri- 
bution does not impact significantly the results. Once all the pa- 
rameters are chosen for a planet, its projected separation on the 
sky plane is calculated taking into account the distance of the 
star, all the orbital parameters and the date of the observations. 

The other important observable needed for comparison to 
the detection limits is the luminosity of the planet. This quan- 
tity is directly related to the planet mass and age, and can be 
estimated using evolutionary models. The predictions from evo- 
lutionary models can vary significantly depending on which 
model is considered. In particular, the cold-start (Marley et al. 
2007; Fortney et al. 2008) models predict much fainter plan- 
ets at young ages for a given mass than the hot-start models 
(Chabrier et al. 2000; Baraffe et al. 2003; Burrows et al. 2003). 
Hot-start models generally use arbitrary initial conditions, in- 
dependent of the outcome of the formation phase, assuming 
that they will be "forgotten" after a few millions of years 
(Baraffe et al. 2002). In contrast, cold-start models have been de- 
veloped to use initial conditions calculated from a core-accretion 
mechanism for the formation of the planets. However, the sim- 



plified treatment of the core-accretion shock results in significant 
uncertainties on the predictions of these models (Fortney et al. 
2008). Existing observations are more consistent with predic- 
tions from the hot-start models rather than the cold-start mod- 
els (see Janson et al. 2011, Fig. 4), although the much fainter 
magnitudes of cold-start planets would not have been detectable 
by most observations. Only hot-start models are considered for 
this study, since the brightness levels of cold-start planets are too 
faint for the typical sensitivity of current observations. 

More precisely, the predicted magnitude from evolutionary 
models (Chabrier et al. 2000; Baraffe et al. 2003) - calculated 
in the CH4s and K' filters of NIRI or in the Ks filter of NaCo 
(I. Baraffe and F. Allard, private communication) - are used to 
estimate the planet magnitudes as a function of mass for the 
age of each target star. Only two extensive evolutionary model 
grids are available - DUSTY and COND (Chabrier et al. 2000; 
Baraffe et al. 2003). Of the two evolutionary models, only the 
COND models extend to sufficiently faint magnitudes and low 
planet masses to reach the detection limits of our observations, 
so the analysis in this paper uses the COND models. As an ex- 
ample, the DUSTY and COND mass-magnitude relations at the 
age of 30 Myr are shown in Fig. 6 which illustrates the limita- 
tion of the DUSTY models compared to the sensitivity limit of 
the observations. Figure 6 also plots the estimates of the masses 
of the outer three HR 8799 planets based on dynamical stability 
limits and evolutionary models (Marois et al. 2010b) and com- 
pares the masses to the predictions of the DUSTY and COND 
models. Within the current capacity to perform empirical com- 
parisons, the //-band and Zf-band magnitudes are consistent with 
the COND models, as the range of estimated masses for the 
HR 8799 planets encompasses both the predicted CH4s and K- 
band magnitudes for that mass range. Of the 42 stars in the sam- 
ple, 9 were observed with the CH4s filter and 33 were observed 
with K-band filter. Spectroscopic and photometric observations 
of the HR 8799 planets have measured differences from the the- 
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oretical DUSTY and COND atmospheres and highlighted the 
need for updated models incorporating non-equilibrium chem- 
istry and an improved treatment of clouds (Jansonet al. 2010; 
Bowler et al. 2010b; Barman et al. 201 la; Galicher et al. 201 1), 
however grids of models including these effects are not currently 
available. 

In the MC simulations, the masses of the generated planets 
are converted to absolute magnitudes in the filter matching the 
observations using the evolutionary models, and the signal of 
each planet is then compared to the 5cr detection limit of the 
observations at the appropriate angular separation. Since we are 
using the ID detection limits, the position angle of the simulated 
planets is simply ignored. A planet is considered as detected if 
its observed signal lies above the detection limit at its projected 
position. The completeness is then calculated as the fraction of 
detected planets over the number of generated planets. 

The two new binaries discovered in the course of our sur- 
vey, HIP 42334 and HIP 104365, were included into the anal- 
ysis taking into account the possible disruption induced by the 
stellar binarity. The limiting values for the semimajor axis of 
the planets were computed using the analytical expressions of 
Holman & Wiegert (1999). They define the limiting values that 
the semimajor axis of a planet can reach and still maintain its 
orbital stability, as a function of the mass-ratio and orbital ele- 
ments of the binary. For HIP 42334 and HIP 104365, the ranges 
of semimajor axes where long-term stability is unlikely are re- 
spectively 85-720 AU and 198-1750 AU, so no planets were 
generated within these intervals. 

Finally, for the 13 targets where candidates were detected 
outside of 320 AU but not followed-up by subsequent observa- 
tions, planets generated in the MC simulation with a projected 
separation larger than the closest candidate detected were con- 
sidered as non-detectable. 

5.4. Mean probability of detection 

The overall sensitivity of our observations can be estimated in 
the mass interval 1-14 Mj up and the semimajor axis interval 5- 
320 AU. For this analysis, we perform MC simulations for each 
target assuming no particular distribution for the mass and semi- 
major axis of the simulated planets. At each point on a linear 
grid of semimajor axis and mass, 10 4 planets are simulated fol- 
lowing the procedure described in Sect. 5.3, and their detectabil- 
ity is tested with the detection limits to calculate a completeness 
map for each target of the sample. The mean survey detection 
thresholds are calculated to illustrate the sensitivity levels, but 
the completeness of each target is used in the statistical analysis 
rather than an aggregate sample sensitivity. 

The mean detection probability of the full survey is estimated 
by averaging all the completeness maps for each of the 42 tar- 
gets. Only the deepest observation of each target was used in the 
MC simulations, since the second epoch observations were gen- 
erally only deep enough to redetect the candidates. For the par- 
ticular case of HR 8799, which is the only target with five very 
deep observations, the completeness maps for the five epochs 
were averaged together before being included with the maps of 
the other targets to avoid giving disproportionate weight to that 
specific target. 

The mean probability of detection for the full survey is rep- 
resented in Fig. 7 with contour plots as a function of mass (1- 
30 Mj up ) and semimajor axis (10-2500 AU), and cuts at specific 
masses. For a planet of given mass and semimajor axis, Fig. 7 
describes the probability that we would have detected it from 
our survey data. Observing bright A-stars necessarily results in 



a greater minimum-mass detection limit, due to the increased 
contrast compared to solar-type stars. This is clearly visible on 
Fig. 7, which shows that our data is most sensitive above 10 Mj up 
and at separations of 50 AU or more, reaching the highest prob- 
abilities of detection between 75 and 300 AU. This is consistent 
with the median distance of our sample (50 pc), a saturation ra- 
dius of 0.3-0.5" for the data, and the fact that the sensitivity 
limits usually reach a plateau at separations larger than 3-4". 
The peak sensitivity reaches 45%, 80% and 94% respectively 
for 5, 10 and 14 Mj up . The overall sensitivity to planets similar 
to HR 8799 bed and f3 Pic is low, with a detection probability 
ranging from 1% for /3 Pic b to 30% for HR 8799 b, due to the 
small number of targets covering the appropriate range of orbital 
separations. The capacity to detect planets like /3 Pic b is partic- 
ularly low, and we consider this complication in the frequency 
calculation in Sect. 5.5. As expected, the sensitivity of the sur- 
vey to BDs is much higher, with a detection probability above 
50% between 75 and 300 AU. 

5.5. Estimation of the planet and brown dwarf systems 
frequency 

In this section we publish the first estimation of the frequency 
of exoplanetary systems in the mass interval 3-14 Mj up and the 
semimajor axis interval 5-320 AU for A-stars. For each target in 
j - 1 . . . 42, the individual target completeness map is used to 
estimate the mean value of pj over the 3-14 Mj up and 5-320 AU 
intervals. The values of pj for each target is then used in Eq. 1 to 
calculate the likelihood for values of / between and 1 . Using 
numerial integration of the likelihood with respect to /, the pos- 
terior distribution p(f\{dj}) is derived using Eq. 2. Finally, 95% 
and 68% confidence intervals are calculated using numerical in- 
tegrations of Eq. 4 and 5. 

Figure 8 presents the probability density function of / calcu- 
lated from our observations, using a linear- flat prior, compared 
to a binomial distribution that represent an ideal case in which 
the observations would be sensitive to all planets in the intervals 
5-320 AU and 3-14 M Jup . We have seen from Sect. 5.4 that our 
sensitivity to low-mass, close-in planets is limited, so the cal- 
culated probability density departs significantly from the ideal 
binomial distribution. When considering both detections around 
HR 8799 and /3 Pic, we estimate the planetary systems frequency 
around A-stars to 8. 7^' % at a confidence level of 68% (lcr). 
This value relies on the assumption made in the Monte-Carlo 
simulations (see Sect. 5.3 and 5.4) that the mass distributions of 
the simulated planets is fiat. We investigate constraints on other 
possible distributions in Sect. 5.6. Although this probability is 
relatively high, it is a direct consequence of our decreased sen- 
sitivity to short-period planets with low mass. Upcoming large 
surveys with increased sensitivity to planets in ~5-50 AU orbits 
will provide a more accurate estimation of the frequency over 
the full 5-320 AU range. 

We note that the inclusion of /? Pic in the sample must be con- 
sidered with care. Although it was included as part of the previ- 
ous deep imaging survey by Chauvin et al. (2010), it is important 
to note that their observations were not sensitive enough to detect 
the planet. For that particular target, we reanalyzed the data pre- 
sented by Bonnefoy et al. (2011), which is more sensitive than 
the data of Chauvin et al. (2010). This could potentially cause 
a bias, because it is unknown whether there are similarly sen- 
sitive unpublished datasets for the other A-stars included from 
the Chauvin et al. (2010) survey. An additional bias could oc- 
cur because /3 Pic may have been selected among all potential 
A-stars targets due to indicators correlated with the presence of 
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Fig. 7. Mean probability of detection of our survey as a function of semimajor axis and planet mass (left) and for specific masses, 
given above or below each curve in Mj up , as a function of semimajor axis (right). The mean probability of detection is obtained 
by averaging the completeness maps of all targets of the survey. /3 Pic b (orange square) and the 4 planets of the HR 8799 system 
(red circles) have been overlaid for reference using the current best estimates of their masse and separation. Note that the HR 8799 
planets are represented at their projected physical separation because the true physical separation is not yet known precisely. 



a planet, such as the known warp in the disk. Because of these 
complications, we also present an estimation of the frequency 
considering /3 Pic as a non-detection 3 . In that case, we estimate 
the frequency to 4.3^ ' % at 68% confidence level (4.3^ 7 3 7 % at 
95% confidence level). 

Finally, since our data is very sensitive to the brown dwarf 
regime, we can make an estimation of the brown dwarf systems 
frequency around A-stars including the detection of HR 7329 B 
in the survey. In the intervals of 14-75 Mj up and 5-320 AU, we 
estimate the frequency to 2.8^ " % at a confidence level of 68%. 

5.6. Constraints on giant planets population distributions 

Another interesting aspect of our sample is the possibility to 
place some initial constraints on the distributions of mass and 
semimajor axis for the giant planet population. The approach 
taken by previous studies was to try extrapolating the distribu- 
tions coming from RV surveys, which are usually parameterized 
with power laws. One of the most complete RV study on this 
topic for FGKM stars is the one performed by Cumming et al. 
(2008) using the sensitivity limits from the Keck Planet Search. 
They determine that the planet frequency for solar-types stars is 
10.5% for 0.3-10 Mj up planets with periods less than 1826 days, 
and that the mass and semimajor axis distributions are best mod- 
eled with power laws respectively of index 4 a = —1.31 and 
P = -0.61, i.e. dN oc M~ l3l dM and dN oc a -°- 6l da. For early- 
type stars, Johnson et al. (2010a) find a positive correlation be- 
tween planet frequency and primary mass, and determine a fre- 
quency of 1 1 + 2% for stars in the range 1.3—1 .9 M (see Table 6 
for their measured frequency in different bins of stellar mass). 
For the distributions of mass and semimajor axis around A-stars, 



3 The contrast curve of Chauvin et al. (2010) is used for that estima- 
tion. 

4 We define the power law indices with respect to linear bins rather 
than the logarithmic bins of Cumming et al. (2008). Also, our f} defines 
the semimajor axis distribution rather than the period distribution. 



[3,14] M Jup - [5,320] AU 
i i 1 1 1 i i 



2.5 - C(42,2) 2 detections: -» HR 8799 




Frequency f (%) 

Fig. 8. Representation of the probability density of / given the 
observations \dj) (posterior distribution) as a function of the 
planetary systems frequency /, in the semimajor axis interval 5- 
320 AU and mass interval 3-14 Mj up , for two detections around 
HR 8799 and fi Pic (see Sect. 5.5 for details). The posterior dis- 
tribution is obtained from Eq. 2 assuming a linear-flat prior for 
the frequency and a flat mass distribution for the planets. An 
equal-tail CI has been calculated using Eq. 4 and 5 for confi- 
dence levels a of 68% (light grey area) and 95% (dark grey area). 
The dashed lines gives the binomial distributions for one or two 
detections out of 42 targets, which represent the probability den- 
sity that would be obtained in an ideal case where the observa- 
tions would be sensitive to all planets in the range 5-320 AU and 
3-14 M, up . 
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Fig. 9. Contour plots showing the confidence with which we can 
exclude a model with specific parameters for the mass distribu- 
tion power law (a, from -1.5 to 1.5, top to bottom), the semima- 
jor axis distribution power law index (/?, from -2.5 to 0.2) and 
the semimajor axis distribution power law cutoff (a C utoff, from 
10 to 320 AU). The hatched grey area corresponds to the values 
of flcutoft within the separation of the known planet HR 8799 b 
(dashed line). Given the existence of the planet at a separation 
of ~68 AU, values of a cuto ff below that separation are de facto 
excluded. The best fit values of Cumming et al. (2008) for solar- 
type stars in terms of mass and semimajor axis are a = —1.31 
and /3 = -0.61. Bowler et al. (2010a) report a 50% confidence 
for values of j3 higher than zero. 



Bowler et al. (2010a) use the detections of their survey to place 
constraints on the possible values of a and /}. They find that the 
values of Cumming et al. (2008) are inconsistent with their ob- 
servations, with a probability below 1 %, and consequently they 
favor positive values for both parameters. 

For our study we assume the frequency / is known over a 
certain range of mass and semi-major axis from RV surveys, 
and we model the mass and semimajor axis distributions as 
power laws. A cutoff on the semimajor axis distribution (a C utoff) 
is added to define an upper limit to the separation at which plan- 
ets are present. To obtain meaningful results it is important to 
define the range of mass and semimajor axis over which the 
frequency / is known and valid. In our case, we use the value 
reported by Bowler et al. (2010a) and Johnson et al. (2010a) in 
their survey of giant planets around old A stars: / = 11+2% 
measured for planets in the ranges 0.5-14 Mj up and 0.1-3.0 AU 
(median sensitivity of their RV data). 

However, our observations are not sensitive to the same 
range of mass and semimajor axis as the RV data, so the value 
of / needs to be scaled to be valid inside the range of parame- 
ters where we are sensitive, [m^, m max ] for the planet mass and 
[a m i n , a ma x] for the semimajor axis. Once the indices a and f3 are 
chosen, the power laws can be integrated over [0.5,14] Mj up and 
[0.1,3.0] AU for the RV, and over [m m i n ,m max ] and [a m in,amax] 
for the direct imaging. Finally, the value of / is scaled by the ra- 
tio of the direct imaging integrated values over the RV integrated 
values. This normalization ensures that the new frequency, /', is 
valid over [m m i a ,m mia ] and [amin,flmax], and that its value over 
[0.5,14] M Jup and [0.1,3.0] AU is always equal to / = 11% to 
match the RV data. In our case, we used the current knowledge 
of the properties of the detected companions around HR 8799 
and /3 Pic to define the appropriate range of mass and semima- 
jor axis: intervals of [3.5,12] Mj up and [8,68] AU were used to 
encompass these currently known detections. 

For each of the targets in our sample, we simulate different 
populations of planets with distributions generated on a grid of 
values for a (-1.5 to 1.5 by steps of 0.1),/? (-2.5 to 0.5 by steps 
of 0.1) and the semimajor axis distribution cutoff, a cutc .ff (10 to 
320 AU by steps of 10 AU). For each point of the grid, 10 5 plan- 
ets are simulated and the fraction of detectable planets in the 
intervals [3.5,12] Mj up and [8,68] AU is measured. If we assume 
that all stars of the sample have one planet in that range, the sum 
of these fractions over our 42 stars would give the number of de- 
tection that would be expected from the observations (taking into 
account the sensitivity around each target). However, we have as- 
sumed that the frequency of planetary systems is equal to /', so 
the sum needs to be multiplied by /' to obtain the expected num- 
ber of detections over [3.5,12] Mj up and [8,68] AU, taking into 
account the known frequency from RV surveys. Finally, we can 
use Poisson statistics to obtain the probability that this expected 
number of detections matches our real detections. In Fig. 9 we 
represent contour plots showing the confidence with which we 
can exclude different combinations of parameters given the de- 
tections in the range [3.5,12] Mj up and [8,68] AU. Before de- 
scribing the results, we remind that increasing the value of a 
is equivalent to increasing the proportion of high-mass planets, 
and that increasing the value of /3 is equivalent to increasing the 
proportion of planets at large orbital separations. We first note 
that given the detection of the planet HR 8799 b at an orbital 
separation of 68 AU, we assume that the cutoff for the semi- 
major axis distribution, a cu toff, lies at higher orbital separations. 
In the three contour plots, the combination of parameters that 
can be excluded with 95% confidence defines a vertical area for 
flcutoff ^ 68 AU. The general trend when increasing a, i.e. when 
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Fig. 10. Comparison of the mean probability of detection for 
the present survey (plain orange contours) with a simulated sur- 
vey with GPI of the same 42 targets (dashed blue contours). 
The mean probability of detection is obtained using Monte- 
Carlo simulations as described in Sect. 5.4. In addition to 
/? Pic b (orange square) and the HR 8799 planets (red circles), 
the giant planets detected around old A stars by Johnson et al. 
(2010b, 201 1) have been overplotted (purple triangles). The dot- 
ted purple line shows the median detection threshold of the 
Bowler et al. (2010a) RV survey around old A-type stars. The 
HR 8799 planets are represented at their projected physical sep- 
aration because the true physical separation is not yet known 
precisely. Note that Fomalhaut b has not been included due to its 
uncertain nature (see e.g. Janson et al. 2012). 

increasing the proportion of high-mass planets, is to move the 
allowed range of parameters towards smaller values for /3. 

The top plot of Fig. 9 is close to the value of a = -1.31 
reported by Cumming et al. (2008). With this value of a, their 
value of p = -0.61 is not in the range of parameters allowed 
by our observations. When exploring higher values of a, as sug- 
gested by Bowler et al. (2010a), a value of -0.61 for (3 becomes 
even less likely given a cutoff of fl C utoff ^ 68 AU. We also note 
that Bowler et al. (2010a) report higher confidence for positive 
values of both a and However, when increasing a in our sim- 
ulations, the positive values of /? become less likely, suggesting 
a possible change in the population of exoplanets between the 
short separation planets detected by RV and the planets at very 
wide orbital separations around A-stars. 

6. Discussion and conlusions 

The present survey results on A-star wide orbit planetary sys- 
tems frequency and population distribution can be compared 
and combined with previous results to investigate the effects of 
host star mass, the predictions of theoretical models, the impact 
of debris disks, and the current state of the planet population 
across a range of orbital separations. Table 6 summarizes the 
frequencies and survey sensitivities reported in both AO (cur- 
rent work and Nielsen & Close 2010, which combines the results 
of Masciadri et al. 2005, Biller et al. 2007, and Lafreniere et al. 
2007a) and RV surveys (Johnson et al. 2010b). These studies 
were selected for comparison, since each covers a large sample 



and considers the same 68% confidence interval. For the sam- 
ple with the most similar target mass range, the evolved A-stars, 
the AO measurement of the wide orbit planetary systems fre- 
quency range of 5.9-18.8% encompasses the close orbit planet 
frequency of 1 1 ± 2%. Direct imaging searches covering similar 
orbital separations to this study, but targeting lower mass stars 
(Nielsen & Close 2010) find an upper limit of <20%, while RV 
searches for close orbit planets measure a declining frequency 
for the successively lower mass samples (Johnson et al. 2010a). 
From the current statistics on the wide orbit imaged planet pop- 
ulation, it is not possible to determine how the planetary sys- 
tems frequency scales with host star mass. Theoretical mod- 
els have predicted a rising (Kennedy & Kenyon 2008), peaked 
(Ida & Lin 2005), or declining (Kornet et al. 2006) planet fre- 
quency on mass, with models incorporating different treatments 
of factors such as the location and evolution of the snow line, 
the initial disk size and its dependence on host star mass, and the 
orbital migration of planets. 

Since the survey was designed to detect planets, the full sam- 
ple is extremely sensitive to brighter brown dwarf companions 
in wide orbits (probability of detection > 90% in -75-300 AU, 
see Fig. 7), enabling an initial assessment of the brown dwarf 
systems frequency around stars more massive than the Sun. 
Among the 42 targets, one brown dwarf was previously iden- 
tified around HR 7329 (Lowrance et al. 2000; Neuhauser et al. 
2011), yielding a brown dwarf systems frequency of 2.8^^% 
(see Sect. 5.5). This low level of brown dwarf companions is 
consistent with the 3.2 3 2 1 7 % frequency measured for a larger 
scale survey of young (3-3000 Myr), solar-type (F5-K5) stars 
covering a wider range of separations (28-1590 AU) with an 
AO survey (Metchev & Hillenbrand 2009). 

Theoretical models of the formation of giant planets by grav- 
itational instability have proposed that the directly imaged plan- 
ets may be the lowest mass products of the same process that 
forms binary stars and brown dwarfs (Kratteret al. 2010), and 
these simulations predict a higher proportion of brown dwarfs 
than planets. While the statistics of the current study are limited, 
the similarity of the brown dwarf and planetary systems frequen- 
cies suggests that the imaged planets may have formed through a 
different process such as core accretion (e.g. Pollack et al. 1996). 
An investigation into the conditions for gravitational instability 
for the specific case of the HR 8799 planets (Nero & Bjorkman 
2009) found that the likelihood of this formation mechanism in- 
creased with radius towards the outermost planets. Other pos- 
sible explanations for the wide orbits of some of the directly 
imaged planets include scattering to larger orbital radii from 
interactions in multiple planet systems (e.g. Veras & Armitage 
2004), outward migration in a resonance with multiple planets 
(e.g. Crida et al. 2009), and outward migration through a planet- 
disk interaction (e.g. Veras & Armitage 2004). 

The full sample of 42 young, A- and F-stars includes 
17 systems with dusty debris disks sustained by the ongo- 
ing collisional grinding of planetesimals into smaller particles 
(Backman & Paresce 1991) or an event such as a catastrophic 
collision of planets (Cameron 1997; Melis et al. 2010). The 3 
targets in this sample with imaged planets or brown dwarfs all re- 
side in systems encircled by dust disks. Additionally, Fomalhaut 
(not included in this study) would provide another example of a 
planet-disk system (Kalas et al. 2008) if the presence of a planet 
was confirmed (see e.g. Janson et al. 2012). While not all A-stars 
with debris disks harbor massive giant planets in wide orbits, the 
frequency of planetary systems appears higher among the targets 
with excess emission from dust. The combination of resolved 
disk maps and planet imaging is an especially powerful tool to 
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Table 6. Samples and methods comparison. 



Sample 


Host mass 


Tech. 


Frequency 


Sep. range 


Planet mass limit 


Reference 




(M Q ) 






(AU) 


(M Jup ) 




F5-A0 


-1.5-3.0 


AO 


5.9-18.8% 


5-320 


3-14° 


this work 


Evolved A 


1.3-1.9 


RV 


1 1 ± 2% 


0.1-3 


>0.2-1.3* 


Johnson et al. (2010a) 


K7-F2 


0.7-1.5 


AO 


< 20% 


25-856 


>4 


Nielsen & Close (2010) 


FGK 


0.7-1.3 


RV 


6.5 ± 0.7% 


0.01-3 


>0.5-0.9 c 


Johnson et al. (2010a) 


M5-M0 


0.2-0.6 


AO 


< 20% 


9-207 


>4 


Nielsen & Close (2010) 


M 


0.1-0.7 


RV 


2.5 ± 0.9% 


0.01-3 


>0. 1-0.5'' 


Johnson et al. (2010a) 



Notes. {a) The detection limit is a function of separation, as detailed in Fig. 7. The frequency estimation is based on a flat distribution for the 
planets mass. < * ) The detection limit is a function of separation, as reported in Bowler et al. (2010a). (c) Scaled from the limits in Bowler et al. 
(2010a) using the fixed RV amplitude cutoff of K > 20 m s _1 reported in Johnson et al. (2010a) and the ratio of the average mass of the bin to the 
average mass of the evolved A-stars, taken to the power of 2/3. 



understand the planet-disk interactions that may sculpt young 
planetary systems (e.g. Liou & Zook 1999; Kuchner & Holman 
2003; Wyatt 2006; Quillen & Faber 2006). The structure of the 
Fomalhaut disk inner edge has been compared with dynamical 
models of planets with different masses (Chiang et al. 2009), and 
the size and shape of the HR 8799 disk have been compared 
with models of the dynamically cleared zones and orbital mi- 
gration history (Patience et al. 2011). Advances in disk imag- 
ing with ALMA and planet imaging with extreme AO will en- 
able more detailed studies to investigate planet-disk interactions 
which may generate structures such as asymmetries and clumps. 

By combining the results and sensitivities of existing di- 
rect imaging A-star planet searches with those of RV pro- 
grams targeting retired A-stars, a comprehensive summary of 
the currently known A-star planet population is given in Fig. 10. 
Analysis of the current data already shows distinct differences. 
The close orbit A-star planets are best fit by a distribution that 
is flat or rising with increasing orbital separation (Bowler et al. 
2010a) and rising with increasing planet mass (Bowler et al. 
2010a; Johnson et al. 2010a). In contrast, the wide orbit planet 
data, incorporating the outer cutoff implied by the outermost 
HR 8799 planet, are consistent with a distribution that is declin- 
ing with increasing orbital separation. This result is independent 
of whether the distribution of planets is rising, flat, or declining 
with planet mass, as shown in the three panels of Fig. 9 (the al- 
lowed region always has a negative power law index). Thus, the 
existing surveys have identified the boundaries of planets around 
A-stars - from -0.6 AU (Bowler et al. 2010a) to -70 AU - and 
have indicated that the data cannot be fit by a single distribution 
of properties. 

Figure 10 also presents the sensitivity of a simulated 
survey with the upcoming high-contrast imager Gemini/GPI 
(Macintosh et al. 2008). Using simulated contrast curves for that 
instrument, we performed Monte-Carlo simulations similar to 
the ones described in Sect. 5.4 on the 42 targets of our sample. 
They show that extreme AO imaging will provide the crucial 
coverage to connect the close orbit and wide orbit populations 
and to reveal the full distribution of planets as a function of sep- 
aration and planet mass. Extreme AO data will determine if the 
transition from a flat or rising population of close orbit planets to 
a declining population of wide orbit planets is smooth or discon- 
tinuous and may discover a peak in the separation distribution. 
The total number of planets detected in extreme AO surveys, 
combined with trends measured for close orbit A-star planets, 
can be used as a test of formation models (e.g. Crepp & Johnson 
2011). Since the upcoming extreme AO instrument include in- 
tegral field units, it will be possible to investigate the nature of 
the exoplanet atmospheres in addition to the population statis- 



tics. Spectra of the currently imaged planetary mass companions 
have revealed differences with brown dwarfs of similar tempera- 
tures (e.g. Janson et al. 2010; Patience et al. 2010; Barman et al. 
2011a,b; Skemeretal. 2011), and upcoming AO observations 
will further explore the architectures and atmospheres of exo- 
planets around A-stars. 
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Appendix A: Point sources properties 

Table A.l summarizes the properties of all point sources de- 
tected with our observations. 
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Table A.l. Properties of all detected candidate companions at the different epochs. 



Name 


CC# 


Date 


Epoch 


Instrument 


Filter 


Sep. 


PA 




Amag 


Proj. sep. 


Status" 








(yr) 








(as) 


(deg) 


(mag) 


(AU) 




HIP 104365 





2007-09-18 


2007.71 


NIRI 


CH4s 


9.545 ±0.021 


15.27 ± 


0.13 


7.4 


± 0.1 


526.2 ± 1.2 


Comoving 


HIP 104365 





2008-09-09 


2008.69 


NIRI 


CH4s 


9.543 ±0.021 


15.25 ± 


0.13 


7.4 


±0.1 


526.1 ± 1.2 


Comoving 


HIP 104365 


1 


2007-09-18 


2007.71 


NIRI 


CH4s 


9.595 ±0.021 


14.27 ± 


0.13 


7.5 


±0.1 


529.0 ± 1.2 


Comoving 


HIP 104365 


1 


2008-09-09 


2008.69 


NIRI 


CH4s 


9.601 ±0.021 


14.20 ± 


0.13 


7.2 


±0.1 


529.3 ± 1.2 


Comoving 


HIP 104365 


2 


2007-09-18 


2007.71 


NIRI 


CH4s 


13.131 ±0.021 


62.07 ± 


0.09 


14.9 


±0.1 


723.9 ± 1.2 


Undefined 


HIP 10670 





2007-09-16 


2007.71 


NIRI 


CH4s 


4.650 ±0.006 


253.27 ± 


0.07 


15.6 


±0.1 


160.1 ±0.2 


Background 


HIP 10670 





2008-10-12 


2008.78 


NIRI 


CH4s 


4.665 ± 0.006 


253.93 ± 


0.07 


14.4 


±0.2 


160.7 ±0.2 


Background 


HIP 10670 


1 


2007-09-16 


2007.71 


NIRI 


CH4s 


10.668 ±0.021 


161.63 ± 


0.11 


13.7 


±0.1 


367.4 ±0.7 


Background 


HIP 10670 


1 


2008-10-12 


2008.78 


NIRI 


CH4s 


10.599 ±0.021 


161.88 ± 


0.12 


12.7 


±0.1 


365.0 ±0.7 


Background 


HIP 10670 


2 


2007-09-16 


2007.71 


NIRI 


CH4s 


13.844 ±0.021 


344.29 ± 


0.09 


13.5 


±0.1 


476.8 ±0.7 


Ambiguous 


HIP 10670 


2 


2008-10-12 


2008.78 


NIRI 


CH4s 


13.982 ±0.021 


344.42 ± 


0.09 


12.4 


±0.1 


481.5 ±0.7 


Ambiguous 


HIP 110935 





2009-06-24 


2009.48 


NaCo 


A s 


0.582 ±0.007 


52.74 ± 


0.66 


9.6 


±0.7 


25.7 ±0.3 


Undefined 


HIP 11360 





2007-12-30 


2007.99 


NIRI 


CH4s 


12.628 ±0.021 


194. 10 ± 


0.10 


14.5 


±0.1 


571.2 ±0.9 


Background 


HIP 11360 





2008-10-12 


2008.78 


NIRI 


CH4s 


12.572 ±0.021 


194.76 ± 


0.10 


12.8 


±0.3 


568.6 ±0.9 


Background 


HIP 14551 





2009-12-21 


2009.97 


NaCo 


A s 


6.058 ±0.007 


28.65 ± 


0.06 


14.9 


±0.2 


331.0±0.4 


Undefined 


HIP 15648 





2008-10-17 


2008.79 


NIRI 


K' 


11.528 ±0.021 


225.35 ± 


0.11 


10.4 


±0.1 


532.2 ± 1.0 


Undefined 


HIP 16449 





2007-09-22 


2007.72 


NIRI 


CH4s 


10.193±0.021 


25.77 ± 


0.12 


11.5 


±0.1 


733.3 ±1.5 


Background 


HIP 16449 





2008-10-12 


2008.78 


NIRI 


CH4s 


10.184±0.021 


25.46 ± 


0.12 


9.2 


±0.1 


732.6 ±1.5 


Background 


HIP 22226 





2008-01-15 


2008.04 


NIRI 


CH4s 


12.156±0.021 


148.90 ± 


0.10 


14.7 


±0.5 


975.6 ±1.7 


Background 


HIP 22226 





2008-11-18 


2008.88 


NIRI 


CH4s 


12. 114 ±0.021 


148.97 ± 


0.10 


12.4 


±0.3 


972.3 ± 1.7 


Background 


HIP 23296 





2010-01-03 


2010.01 


NIRI 


A" 


6.510 ±0.021 


308.63 ± 


0.19 


13.8 


±0.1 


322.9 ± 1.0 


Undefined 


HIP 23296 


1 


2010-01-03 


2010.01 


NIRI 


A" 


8.016 ±0.021 


99.07 ± 


0.15 


15.1 


±0.1 


397.6 ±1.0 


Undefined 


HIP 23296 


2 


2010-01-03 


2010.01 


NIRI 


K' 


9.174 ±0.021 


329.02 ± 


0.13 


7.5 


±0.1 


455.0 ±1.0 


Undefined 


HIP 23296 


3 


2010-01-03 


2010.01 


NIRI 


K' 


9.727 ±0.021 


29.27 ± 


0.13 


15.2 


±0.1 


482.5 ± 1.0 


Undefined 


HIP 23296 


4 


2010-01-03 


2010.01 


NIRI 


K' 


11.360 ±0.021 


173.73 ± 


0.11 


13.0 


±0.1 


563.5 ±1.0 


Undefined 


HIP 23296 


5 


2010-01-03 


2010.01 


NIRI 


K' 


13.586 ±0.021 


73.30 ± 


0.09 


12.3 


±0.1 


673.9 ±1.0 


Undefined 


HIP 26624 





2008-11-14 


2008.87 


NIRI 


K' 


8.753 ±0.021 


232.98 ± 


0.14 


11.0 


±0.1 


372.6 ±0.9 


Undefined 


HIP 26624 


1 


2008-11-14 


2008.87 


NIRI 


K' 


11.914 ±0.021 


42.62 ± 


0.10 


13.1 


±0.1 


507.2 ±0.9 


Undefined 


HIP 34782 





2009-12-21 


2009.97 


NaCo 


K s 




8.667 ±0.007 


316.66 ± 


0.04 


14.2 


±0.3 


412.1 ±0.3 


Background 


HIP 34782 





2012-01-11 


2012.03 


NaCo 


K, 




8.544 ± 0.007 


317.54 ± 


0.05 


14.0 


±0.3 


406.3 ±0.3 


Background 


HIP 34782 


1 


2012-01-11 


2012.03 


NaCo 


A s 




8.552 ±0.007 


145.51 ± 


0.05 


12.2 


±0.1 


406.6 ±0.3 


Undefined 


HIP 35567 





2009-12-20 


2009.97 


NaCo 


K f 




4.193 ±0.007 


7.74 ± 


0.09 


10.1 


±0.1 


297.8 ±0.5 


Background 


HIP 35567 





2011-12-18 


2011.96 


NaCo 


a; 




4.097 ± 0.007 


8.35 ± 


0.09 


8.9 


±0.1 


291.0 ±0.5 


Background 


HIP 35567 


1 


2009-12-20 


2009.97 


NaCo 


a; 




4.498 ± 0.007 


61.62 ± 


0.09 


10.2 


±0.1 


319.4 ±0.5 


Background 


HIP 35567 


1 


2011-12-18 


2011.96 


NaCo 


a; 




4.467 ± 0.007 


63.31 ± 


0.09 


9.1 


±0.1 


317.2 ±0.5 


Background 


HIP 35567 


2 


2009-12-20 


2009.97 


NaCo 






5.499 ±0.007 


8.61 ± 


0.07 


8.9 


±0.1 


390.5 ±0.5 


Background 


HIP 35567 


2 


2011-12-18 


2011.96 


NaCo 






5.407 ±0.007 


9.14± 


0.07 


7.5 


±0.1 


384.0 ±0.5 


Background 


HIP 41307 





2009-12-21 


2009.97 


NaCo 


A s 




5.792 ±0.007 


86.15 ± 


0.07 


15.7 


±0.2 


217.3 ±0.3 


Background" 


HIP 42334 





2009-12-20 


2009.97 


NaCo 


A s 




2.010 ±0.007 


80.29 ± 


0.19 


11.8 


±0.1 


142.9 ±0.5 


Background 


HIP 42334 





2011-04-27 


2011.32 


NaCo 


a; 




2.047 ± 0.007 


80.39 ± 


0.19 


12.0 


±0.2 


145.5 ±0.5 


Background 


HIP 42334 


1 


2009-12-20 


2009.97 


NaCo 


A f 




3. 176 ±0.007 


1 1.52 =t 


0.12 


8.2 


±0.1 


225.7 ±0.5 


Comoving 


HIP 42334 


1 


2011-04-27 


2011.32 


NaCo 


A f 




3.166 ±0.007 


11.51 ± 


0.12 


8.3 


±0.1 


225.0 ±0.5 


Comoving 


HIP 42334 


2 


2009-12-20 


2009.97 


NaCo 


a; 




5.291 ±0.007 


288.05 ± 


0.07 


12.9 


±0.1 


376.0 ±0.5 


Background 


HIP 42334 


2 


2011-04-27 


2011.32 


NaCo 


a; 




5.265 ±0.007 


288.37 ± 


0.07 


12.9 


±0.1 


374.2 ±0.5 


Background 


HIP 42334 


3 


2009-12-20 


2009.97 


NaCo 


a; 




5.616 ±0.007 


290.18 ± 


0.07 


15.3 


±0.3 


399.1 ±0.5 


Background 


HIP 42334 


3 


2011-04-27 


2011.32 


NaCo 


A s 


5.584 ±0.007 


290.52 ± 


0.07 


14.9 


±0.6 


396.9 ±0.5 


Background 


HIP 44923 





2008-03-23 


2008.22 


NIRI 


CH4s 


13.860±0.021 


236.78 ± 


0.09 


14.1 


±0.1 


1159.8 ±1.8 


Undefined 



Table A.l. continued. 



Name 


CC# 


Date 


Epoch 


Instrument 


Filter 


Sep. 


PA 




Amag 


Proj. sep. 


Status" 








(yr) 






(as) 


(deg) 


(mag) 


(AU) 






HIP 53771 





2010-03-06 


2010.18 


NaCo 


K s 


5.800 + 0.007 


266.75 ± 


0.07 


13.5 


±0.3 


354.3 ± 


0.4 


Background 


HIP 53771 





2012-01-12 


2012.03 


NaCo 




5.731 ±0.007 


267.64 ± 


0.07 


13.3 


±0.1 


350.1 ± 


0.4 


Background 


HIP 53771 


1 


2010-03-06 


2010.18 


NaCo 


a; 


6. 120 ±0.007 


251.12 ± 


0.06 


11.7 


±0.1 


373.9 ± 


0.4 


Background 


HIP 53771 


1 


2012-01-12 


2012.03 


NaCo 


a s 


6.058 ±0.007 


251.73 ± 


0.06 


11.7 


±0.1 


370.1 ± 


0.4 


Background 




2 


zUlz-Ul-lz 


2U12.U3 


NaCo 


K s 


o.D4j ±U.U(J/ 


144.52 ± 


0.05 


12.5 


±0.2 


522.0 ± 


0.4 


Undefined 


tttt» no^n 

HIP 57328 





2008-03-24 


2008.23 


NIRI 


CH4s 


12.505 ± 0.021 


352.40 ± 


0.10 


16.8 


±0.2 


467.8 ± 


0.8 


Undefined 


HIP 61960 





2008-03-21 


2008.22 


NIRI 


CH4s 


12.732 ±0.021 


57.65 ± 


0.10 


15.6 


±0.2 


461.8 ± 


0.8 


Undefined 


HIP 66634 





2009-02-03 


2009.09 


NIRI 


K' 


6.955 ±0.021 


230.69 ± 


0.18 


14.1 


±0.2 


372.7 ± 


1.1 


Undefined 


HIP 69732 





2009-02-02 


2009.09 


NIRI 


K' 


11.462 ±0.021 


319.19 ± 


0.11 


12.3 


±0.1 


348.0 ± 


0.6 


Undefined 


HIP 78078 





2009-06-24 


2009.48 


NaCo 


A s 


6.670 ±0.007 


76.68 ± 


0.06 


13.3 


±0.2 


340.8 ± 


0.4 


Undefined 


HR 7329 





2008-08-07 


2008.60 


NaCo 


H 


4.185 ±0.007 


161.61 ± 


0.09 


< £ 


.12 


201.8 ± 


0.3 


Comoving 1, 


HIP 99273 





2007-08-10 


2007.61 


NIRI 


CH4s 


4.894 ±0.006 


140.41 ± 


0.07 


15.1 


±0.1 


255.6 ± 


0.3 


Background 


HIP 99273 





2008-09-04 


2008.68 


NIRI 


CH4s 


4.802 ±0.006 


140.26 ± 


0.07 


15.3 


±0.2 


250.8 ± 


0.3 


Background 


HIP 99273 


1 


2007-08-10 


2007.61 


NIRI 


CH4s 


7.577 ±0.021 


74.29 ± 


0.16 


16.4 


±0.3 


395.7 ±1.1 


Background 


HIP 99273 


1 


2008-09-04 


2008.68 


NIRI 


CH4s 


7.566 ±0.021 


73.65 ± 


0.16 


16.4 


±0.3 


395.1 ±1.1 


Background 


HIP 99273 


2 


2008-09-04 


2008.68 


NIRI 


CH4s 


8.089 ±0.021 


352.47 ± 


0.15 


16.1 


±0.3 


422.4 ± 


1.1 


Background 


HIP 99273 


2 


2007-08-10 


2007.61 


NIRI 


CH4s 


8.005 ±0.021 


352.70 ± 


0.15 


16.1 


±0.2 


418.0 ± 


1.1 


Background 


HIP 99273 


3 


2007-08-10 


2007.61 


NIRI 


CH4s 


8.917 ±0.021 


304.99 ± 


0.14 


16.1 


±0.2 


465.6 ± 


1.1 


Background 


HIP 99273 


3 


2008-09-04 


2008.68 


NIRI 


CH4s 


8.986 ±0.021 


305.22 ± 


0.14 


16.1 


±0.3 


469.2 ± 


1.1 


Background 


HIP 99273 


4 


2007-08-10 


2007.61 


NIRI 


CH4s 


10.497 ±0.021 


147.75 ± 


0.12 


13.1 


±0.1 


548.2 ± 


1.1 


Background 


HIP 99273 


4 


2008-09-04 


2008.68 


NIRI 


CH4s 


10.407 ±0.021 


147.77 ± 


0.12 


13.0 


±0.1 


543.5 ± 


1.1 


Background 


HIP 99273 


5 


2007-08-10 


2007.61 


NIRI 


CH4s 


10.841 ±0.021 


220.66 ± 


0.11 


8.8 


±0.1 


566.1 ± 


1.1 


Background 


HIP 99273 


5 


2008-09-04 


2008.68 


NIRI 


CH4s 


10.845 ±0.021 


221.04 ± 


0.11 


8.5 


±0.1 


566.3 ± 


1.1 


Background 


HIP 99273 


6 


2007-08-10 


2007.61 


NIRI 


CH4s 


12.318 ±0.021 


75.04 ± 


0.10 


14.1 


±0.9 


643.2 ± 


1.1 


Background 


HIP 99273 


6 


2008-09-04 


2008.68 


NIRI 


CH4s 


12.294 ±0.021 


74.66 ± 


0.10 


14.5 


±0.1 


642.0 ± 


1.1 


Background 


HR 8799 





2008-09-05 


2008.68 


NIRI 


CH4s 


0.630 ±0.006 


197.39 ± 


0.55 


11.6 


±0.3 


24.8 ± 


0.2 


Comoving' 7 


HR 8799 


1 


2007-10-17 


2007.79 


NIRI 


CH4s 


0.943 ± 0.006 


314.08 ± 


0.37 


13.0 


±0.4 


37.2 ± 


0.2 


Comoving' 7 


HR 8799 


1 


2008-09-01 


2008.67 


NIRI 


CH4s 


0.943 ± 0.006 


315.37 ± 


0.37 


12.4 


±0.3 


37.2 ± 


0.2 


Comoving'' 


HR 8799 


1 


2008-09-05 


2008.68 


NIRI 


CH4s 


0.945 ± 0.006 


315.28 ± 


0.36 


11.1 


±0.1 


37.2 ± 


0.2 


Comoving' 7 


HR 8799 


1 


2008-10-10 


2008.77 


NIRI 


CH4s 


0.937 ± 0.006 


315.55 ± 


0.37 


12.5 


±0.5 


36.9 ± 


0.2 


Comoving'' 


HR 8799 




2008-10-14 


2008.79 


NIRI 


CH4s 


0.943 ± 0.006 


314.82 ± 


0.37 


10.7 


±0.2 


37.2 ± 


0.2 


Comoving 


HR 8799 


2 


2007-10-17 


2007.79 


NIRI 


CH4s 


1.721 ±0.006 


60.70 ± 


0.20 


13.8 


±0.2 


67.8 ± 


0.2 


Comoving' 7 


HR 8799 


2 


2008-09-01 


2008.67 


NIRI 


CH4s 


l.709± 0.006 


61.70± 


0.20 


13.5 


±0.2 


67.3 ± 


0.2 


Comoving 


HR 8799 


2 


2008-09-05 


2008.68 


NIRI 


CH4s 


1.722 ±0.006 


61.11 ± 


0.20 


12.2 


±0.1 


67.8 ± 


0.2 


Comoving'' 


HR 8799 


2 


2008-10-10 


2008.77 


NIRI 


CH4s 


1.768 ±0.006 


60.99 ± 


0.20 


14.1 


±0.4 


69.7 ± 


0.2 


Comoving' 7 


HR 8799 


2 


2008-10-14 


2008.79 


NIRI 


CH4s 


1.731 ±0.006 


59.41 ± 


0.20 


11.8 


±0.1 


68.2 ± 


0.2 


Comoving' 7 


HIP 116431 





2007-09-12 


2007.70 


NIRI 


CH4s 


8.748 ±0.021 


28.89 ± 


0.14 


15.7 


±0.2 


598.8 ± 


1.4 


Background 


HIP 116431 





2008-09-24 


2008.73 


NIRI 


CH4s 


8.761 ±0.021 


28.51 ± 


0.14 


15.6 


±0.3 


599.7 ± 


1.4 


Background 



> 
< 



o 

— 

re 

s 

■5 

o 
— I 

re 
I 



> 



Notes. (fl) Note that in general, only candidates with a projected separation larger than 320 AU were followed-up. " J) Background status confirmed by Janson et al. (2011). <c) Known brown-dwarf 
companion redetected in our data (Lowrance et al. 2000). (<,) Co-motion of the HR 8799 planets has been confirmed by multiple authors (e.g. see Marois et al. 2008b, 2010b; Bowler et al. 2010b; 
Currie et al. 2011). 



